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ABSTRACT 

We present two wide-field (w 5' x 3.5'), diflFraction limited (A/D ~ 0.5" at 10 /^m), broad-band 
10 and 20 /^m images of the Orion Nebula, plus six 7-13 iiva. narrow-band (A/AA ~ 1) images 
of the BN/KL complex taken at the 3.8m UKIRT^ telescope with the MPIA MAX camera. The 
wide field images, centered on the Trapezium and BN/KL regions, are mosaics of 35" x 35" frames 
obtained with standard chopping and nodding technique and reconstructed using a new restoration 
method developed for this project. They show the filamentary structure of the dust emission from the 
walls of the HII region and reveal a new remarkable group of arc-like structures « 1' to the South of 
the Trapezium. The morphology of the Ney- Allen nebula, produced by wind-wind interaction in the 
vicinity of the Trapezium stars, suggests a complex kinematical structure at the center of the Cluster. 
We find indications that one of the most massive members of the cluster, the B0.5V star ^^Ori-D, is 
surrounded by a photoevaporated circumstellar disk. Among the four historic Trapezium OB stars, 
this is the only one without a binary companion, suggesting that stellar multiplicity and the presence 
of massive circumstellar disks may be mutually exclusive. In what concerns the BN/KL complex, we 
find evidence for extended optically thin silicate emission on top of the deep 10 /im absorption feature. 
Assuming a simple two component model, we map with ~ 0"5 spatial resolution the foreground optical 
depth, color temperature and mid-IR luminosity of the embedded sources. We resolve a conspicuous 
point source at the location of the IRc2-A knot, approximately 0"5 north of the deeply embedded HII 
region "I" . We analyze the spectral profile of the 10 /xm silicate absorption feature and find indication 
for grain crystallization in the harsh nebular environment. In the OMC-1 South region, we detect 
several point sources and discuss their association with the mass loss phenomenology observed at 
optical and millimeter wavelengths. Finally, we list the position and photometry of 177 point sources, 
the large majority of which detected for the first time in the mid-IR. Twenty two of them lack a 
counterpart at shorter wavelengths, and are, therefore, candidates for deeply embedded protostars. 
The comparison of photometric data obtained at two different epochs reveals that source variability 
at 10 /ini is present up to w 1 mag level on a time-scale 2 yr. With the possible exception of a 
pair of OB stars, all point sources detected at shorter wavelengths display 10 /im emission well above 
the photospheric level, that we attribute to disks circumstellar emission. The model of Robberto, 
Beckwith, & Panagia (2002, ApJ. 578, 897) provides the simplest explanation for the observed mid-IR 
excess. 

Subject headings: nebulae: Orion Nebula, Infrared, ISM, Data Analysis — Infrared 



1. INTRODUCTION 

The Great Orion Nebula (M42=NGC1976), with its associated cluster of Pre-Main-Sequence stars, represents the 
nearest {d ~ 450 pc) and best studied episode of massive star formation (Lada & Lada 2003). The stellar cluster, 
distributed on a volume of ~ 3 pc radius (Hillenbrand & Hartmann 1998), contains approximately 2500 stars ~ 10^ yr 
old spanning the mass range 45-0.02 Mq (Hillenbrand 1997; Hillenbrand & Carpenter 2000). There is clear evidence 
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of mass segregation, with the most massive OB stars (the Trapezium system) concentrated in the central s=; 0.15 pc, 
where the stehar density approaches 2 x 10^ stars pc^"^. Ultraviolet radiation from the Trapezium stars, mostly the 
06.5V star O^Oii C, has produced the Great Orion Nebula (O'Dell 2001a, b), a blister cavity excavated on the surface 
of the parental molecular cloud OMC-1. The same UV radiation also photo-ablates the accreting circumstellar disks 
of the nearest stars, producing the characteristic photoionized "proplyds" (Laques & Vidal 1979; O'Dell et al. 1993). 
Proplyds provide the most compelling evidence that the star formation process can be dramatically affected by the 
mutual interaction between cluster members. 

The Orion Nebula cluster, especially the core region (Trapezium cluster), has been investigated by several authors 
using data collected at visible (Prosser et al. 1985; Hillenbrand 1997; Robberto et al. 2004) and near-IR wavelengths 
(McCaughrean & Stauffer 1994; Hillenbrand & Carpenter 2000; Lucas & Roche 2000; Luhman et al. 2000; Lada 
et al. 2000). Since circumstellar disks emit in the IR the energy they "passively" absorb from their central stars, 
or "actively" produce by the viscous dissipation processes associated to mass accretion, there is a need for similar 
observations at longer wavelengths, in particular in the mid-IR (A ^ 5 — 20 /im), where ground based telescopes can 
reach s\ib-arcsccond spatial resolution. Mid-IR observations have also the potential of discovering objects more deeply 
embedded within the molecular cloud, otherwise invisible at shorter wavelengths (with the possible exception of hard 
X-rays) and eventually identify sites where protostellar contraction is still ongoing. Moreover, mid-IR data may also 
provide complementary information on the structure of the HII region, through the diffuse emission of the dust grains 
and PAH molecules mixed with the ionized gas, or located in the surrounding photodissociation region. 

The first 3.5 fxm image of the Trapezium cluster taken with a panoramic detector was presented, but only briefly 
discussed, by McCaughrean et al. (1996). A more complete study of the Trapezium cluster at this wavelength has been 
published by Lada et al. (2000) . They found that a fraction 80%-85% of stars are surrounded by circumstellar disks, 
in agreement with earlier studies made at near-IR wavelengths, and identify 78 sources with extremely red colors that 
they classify as candidate prostostars. 

At wavelengths A > 3.5 /xm, the region remains largely unexplored. Maps made with single element bolometers in 
the early days of IR astronomy revealed two prominent sources, the BN/KL complex (Becklin & Neugebauer 1967; 
Kleinman & Low 1967) and the Ney- Allen nebula in the vicinity of the Trapezium (Ney & Allen 1969). Both fields 
have been widely investigated with increasing sensitivity and spatial and spectral resolution. Images at 10 and 20 /xm 
made with panoramic detectors have been published by Gezari et al. (1998); Grenhill et al. (2004); Shuping et al. 
(2004) for BN/KL and by Hayward et al. (1994) for the Ney- Allen nebula. Smith et al. (2004) have recently presented 
an image of the OMC-IS region. 

We have undertaken a large observational program on the Orion Nebula using the mid-IR camera MAX on the UKIRT 
telescope. Our research program consists of two parts: first, to observe at 10 /im with relatively long integration times 
the most prominent proplyds, in order to detect and resolve their structure in the mid-IR. Second, to map the entire 
nebular core, i.e. an area of approximately 5' x 5' covering both the Ney- Allen and BN/KL regions, with shorter 
integration times. In this paper we present the results relative to this second study, complemented by narrow-band 
observations of the BN/KL complex. In Section 2 we present our observations and our data reduction strategy, based 
on a reconstruction method originally developed for this project. In Section 3, we present our main results, the 10 and 
20 ^m wide field images and the photometric catalog of 177 point sources detected at 10 /im. In Section 4, we discuss 
the overall structure of the nebula (Section 4.1), the Ney-AUen region (Section 4.2), the BN-KL complex (Section 4.3), 
the OMC-1 South region (Section 4.4), and finally the main properties of our sample of point sources (Section 4.5). 
In Section 5, we summarize our findings. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Data acquisition 

The wide field maps presented in this paper were obtained combining observations made in November 1998 and 
December 2000 at the 3.8m United Kingdom Infrared Telescope (UKIRT) on Mauna Kea, Hawaii. We used MAX^ 
(Robberto & Herbst 1998), a mid-IR imager based on a Rockwell 128 x 128 Si:As Blocked Impurity Band detector opti- 
mized for high background observations. MAX's optics delivered a plate scale of 0'.'26, corresponding to approximately 
35" X 35" field of view. Images have been obtained through broad N-band and Q-band filters. M. Cohen provided de- 
tailed modeling of our photometric system, including typical Mauna Kea transmission, optics throughput, and detector 
response. For an AOV star, the N-band filter has an effective wavelength Aio = 10.1 /im and AA = 3.13 /xm, whereas 
the Q-band filter has A20 = 19.7 /im, and AA = 2.44 /im. The narrow-band images were obtained on December 18, 
2000, through the standard set of silicate filters centered at 7.7, 8.7, 9.8, 10.2, 11.2 and 12.4 /xm (see Robberto & 
Herbst 1998, for the detailed filter parameters). 

UKIRT is equipped with tip-tilt adaptive correction at the secondary mirror (Hawarden et al. 1999). Exploiting 
this capability, MAX routinely reaches diffraction limited performances at A > 8 /xm with \/D ~ 0'.'5. The guide 
system, however, only allows for chopping throws ^ 30" (a maximum of ~ 20"was in fact suggested for the ~ 5 Hz 
chopping frequency we used). Given the extent of the Orion nebula and the MAX sensitivity, we could not chop onto 
empty sky areas for background subtraction. We experimented with the use of sky frames taken several arcminutes 
away, nodding the telescope at a frequency ~ 0.1 Hz, constrained by the telescope pointing and tip-tilt reference star 
acquisition, but this method produced significant losses of sensitivity and observing efficiency. We therefore opted for 

* the instrument has been decommissioned after December 2000. The 20 /xm mosaic represents the very last set of images taken by 
MAX. 
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the conventional chopping and nodding technique, adopting an observing strategy optimized for image post-processing, 
namely the reconstruction method envisioned in Bertero et al. (2000, and references therein). In Table 1, we list the 
main observing parameters: night of the observation, filter, integration time per single frame, chopping frequency, 
total exposure time on the main beam per each nodding cycle and chopping throw amplitude and direction. 

We scanned the survey region both in declination and right ascension using different chopping and nodding ampli- 
tudes. At each step, the telescope was offset (nodding) by an amount equal to the chopping amplitude, i.e. the new 
pointing position of the telescope on the main beam was coincident with the position imaged in the previous offset 
beam. Directions and amplitudes of chopping and nodding were carefully matched and aligned with the orientation 
of the array columns, or rows. Due to a ~ 1° rotation of the detector with respect to the equatorial coordinates, the 
chopping/nodding direction were not exactly coincident with the NS or EW axis. Therefore, the combined rows or 
columns have a slight drift with respect to the scan direction. The setup process and observing strategy is detailed 
further in Bertero et al. (2000). 

Scans were typically ~ 4 arcmin long, with a small side overlay (~ 10 pixels) to maximize field coverage. Since in 
the two runs the camera was mounted on different focal positions of the telescope, images are rotated by 90° between 
the two runs. The rotation of the optical path allows us to easily individuate ghost images produced by the brightest 
sources, BN in particular. 

Observations were generally performed up to an airmass z = 2 under excellent weather conditions. Absolute flux 
calibration in the N band has been obtained by observing at various airmasses the following standard stars: HR 337 
([N]=-2.06), HR 1040 ([N]=2.55), HR 1713 ([N]=0.06), HR 2061 ([N]=-5.13), HR 4069 ([N]=-0.94), HR 2943 ([N]=- 
0.76), HR 4554 ([N]=2.32) and HR 5340 QN]=-3.17). Our absolute flux calibration assumes for a Lyrae [N]=0 with 
a flux F = 1.116 ± 0.015 x 10~^® Wcm~"/xm~^ at Aio, corresponding to 37.9 Jy. Calibration at 20/zm was based 
on HR 1760 ([Q]=-0.30), HR 2061 ([Q]==5.7), and HR 4069 ([Q]=-1.00). The zero magnitude flux in this case is 
F = 7.887 ± 0.011 X 10^^^ Wcm~^/im~^ at A20, corresponding to 10.2 Jy. The narrow-band images of BN were 
calibrated against HR 1040, HR 1047, HR 1713, HR 4069, HR 2943 and HR 5340, with magnitudes and zero points 
given in Table 2. It is clear that the depth of our images is not uniform, since the total integration time in each scan 
is a function of the adopted pattern (amplitude of the chopping throws and offset between adjacent rows/columns), 
airmass, repeated or discarded frames, etc. Moreover, in the final images, part of the field has been observed on more 
than one night. In general, we estimate that the average integration time in the N-band is of order 200 seconds. 
Concerning the narrow-band images, the total integration time is 103.5 s at all wavelengths except the 8.7 fim. band, 
where is 102.5 s. 

As a final remark, we note that we were able to obtain several scans of the region only at 10 /xm. The 20 jim. data 
were all collected in one night, with only one chopping and nodding amplitude. The final mosaic is therefore shallower 
and potentially more affected by reconstruction artifacts. 

2.2. Data reduction 

The raw images were first median filtered to remove a small number (« 15) of dead pixels and occasional extra-noise 

in the frames containing the saturated BN peak. Cleaned frames were then normalized to the same integration time 
and combined to form chopped and nodded pairs. The final combination in the wide area mosaic required a special 
effort, described in the following section. 

2.2.1. Mosaics 

All chopped and nodded pairs taken on the same night were combined into rows or columns and then into a single 
mosaic. This operation was performed using MOSAIC, an IDL package designed for mid-IR data and publicly available 
on the Internet (Varosi & Gezari 1993). The images were combined by integer pixels, rounding the coordinates of the 
FITS header (accurate to ~ O'.'l) to the nearest value. Figure 1 illustrates the mosaicing process using images directly 
grabbed from the MOSAIC output windows. The mosaics contain for each source two negative counterparts at the 
position of the two offset beams. Hence, we needed to reconstruct a true image. 

The restoration method we have used was specifically developed for this project. The theoretical basis have been 
illustrated in a number of papers (Bertero et al. 1998, 1999, 2000), thus we provide here only a brief outline. Since the 
problem is underconstrained, it admits infinit solutions. In particular, the set of solutions whose chopped and nodded 
images are identically zero constitutes the null space, infinite dimensional, of the chopping and nodding operator. Any 
linear combination of these functions can be added to a reconstructed image without affecting the observed, chopped 
and nodded, image. In the null space it is possible to define a metric and apply minimization processes. Our strategy 
is to search for the minimal positive solution, assuming the celestial signal cannot be negative. The implementation of 
the algorithm is based on the Landweber method (Bertero & Boccacci 1998). It is iterative and therefore produces an 
approximate solution at each step. The number of iterations needs to be optimized to control noise propagation, and 
our criterion has been described in Bertero et al. (2000). When applied to real images, it produces excellent results. In 
some situations, however, it may generate artifacts, i.e. artificial ghost images of the stars (appearing as dark images 
over a bright background or as bright images over a dark background) or discontinuities in the backgroiuid level. The 
location and periodicity of the artifacts are related to the values of image size and chopping amplitude through the 
mathematical structure of the problem (Bertero et al. 1999). For our 10 /um mosaic, they have been minimized by 
replacing at their positions the information provided by the reconstructed images obtained with different chopping 
throws and orientation. Another problem is that the inversion method we have used provides less than ideal results 
when the signal in the image to be inverted is negative. In this case, the number of iterations needed to recover 
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the signal may be larger than the number of iterations corresponding to optimal noise reduction. Also in this case, 
these problems have been greatly mitigated at 10 /um by using images taken with different orientation and chopping 
amplitudes. 

2.2.2. Narrow-band im,ages 

The narrow-band images were registered using the peak of the BN object, then combined and trimmed to the same 
size of 101 X 117 pixels, corresponding to 26'.'8 x Sl'.'O. The astrometric solution was derived assuming for the BN peak 
the following coordinates: RA(2000.0)=05''35'"14!10 and DEC(2000.0)=-05°22'22'.'7, a nominal pixel scale of 0'.'26 
constant across the field, a rotation angle of — 0?8 east of north, as measured on the same night for the wide field area 
project. The BN center was set at the pixel position 64.2, 82.2, and measured using a centroid procedure with 5 pixel 
(I'.'S) radius. 

2.2.3. Photometry 

Point sources have been identified by visual inspection of the original (chopped and nodded) images, and then 
confirmed on the final reconstructed image. A few sources, indicated in Table 3, appear very faint only in the final 
mosaic and don't show any counterparts at shorter wavelengths. It is possible that some are artifacts surviving the 
image cleaning processes described in the previous section, and an independent confirmation is needed to make sure 
they are real. On the other hand, the variations of exposure time and airmass across the image may have rendered 
invisible sources otherwise within reach with our instrument in more ideal conditions. 

Aperture photometry was obtained on individual chop/nod pairs using circular apertures with radii between 2 and 
5 pixels, and applying an aperture correction based on the average encircled energy of the standard stars observed 
during the same night. Typically, all standards were used, with the exception of HR 2061 {a Ori), known to be 
extended in the mid-IR (Rinehard et al. 1998) and easily resolved by MAX. An airmass correction was applied to 
each frame. In several cases, the non-uniformity of the background emission dominates the final photometric accuracy. 
For point sources, an error in the sky estimate affects in a predictable way the results of multiaperture photometry 
with curve of growth correction: an overestimate of the sky will cause the source to become fainter with the radius, 
and vice versa. We used this property to tune the local sky value to obtain photometric results as much as possible 
independent of the aperture. Sometimes, however, the variations of the sky brightness underneath the source or 
the extended PSF profile produced unstable results. In these cases, we adopted a different strategy. Using the 
IDL procedure MIN_CURVE_SURF, we estimated the brightness distribution of the sky below the source using the 
minimum curvature surface fitting a sky annulus surrounding the star. The assumption here is that in the presence 
of strongly variable background a higher order interpolation provides a more reliable estimate of the sky distribution 
under the source. The subtraction of this sky frame from the original one provides stellar images that sit on a relatively 
flat background, and therefore, can be measured with greater accuracy. Results obtained with this method turned out 
to depend on the choice of the inner and outer radii of the annulus, typically 5 and 8 pixels. Our final results were 
then obtained by averaging the results provided by a few different, but reasonable on the basis of visual inspection, 
combinations of inner and outer radius. Their scatter has also been added in quadrature to the total uncertainty. 
In general, photometric measurements obtained with different methods coincide within their errors, especially for 
the smaller apertures (1" — 1.5") and within the same run. The imcertainties we quote combine quadratically the 
uncertainties of the zero points (0.12 mag at 7.7 /im, 0.06 mag at 8.7 /im, 0.07 mag at 9.6 /im, 0.06 mag at 10.2 fim, 
0.05 mag at 11.6 /zm, 0.05 mag at 12.4 nm, 0.06 mag in the broad N-band, and 0.11 mag in the Q-band), the source 
photon statistics and sky uniformity, and the scatter between the best estimates obtained with different methods, or 
between different choices of radius/annulus, as explained. 

All sources that had photometric measurements different between the two runs by an amount comparable to, or 
larger than, the quoted 1-cr error were double checked. A direct comparison with nearby sources provides in several 
cases strong evidence for source variability. For sources isolated, or superimposed on a large and complex background 
(observed with different chopping throws/orientation in the two runs), it is not possible to safely discriminate between 
variability and unaccounted photometric errors. Whereas the main photometric catalog has been obtained analyzing 
the single individual images before mosaicing and reconstruction, we have also examined the final reconstructed 
mosaics, both as a sanity check and to find possible evidence of faint sources that may appear only after all frames 
have been stacked together. 

Concerning the overall sensitivity, the final signal-to-noise ratios of our reconstructed mosaics depend on a number 
of factors such as the integration time of each exposure, the number of overlaying exposures contributing at each point, 
airmass, background brightness, chopping direction and noise propagation in the reconstruction algorithm. Therefore, 
the sensitivity and completeness limits vary across the field in a complex way. To provide an indicative value of 
our sensitivity limit, we observe that our catalog of 10 micron sources contains objects down to [A^] ~ 8 mag (see 
Section 4.5). Since all stars brighter that [N] ~ 6.0 mag are well detected in all partial mosaics, we can assume 6.0 mag 
as our completeness limit at 10 /um, at least in the central 2' x 2' of our field. However, according to this criteria, this 
value raises to [N] ~ 5 in the immediate vicinity of the Trapezium, where the nebular emission becomes stronger. This 
because six sources in the range [N] = 5.0 — 6.0 have been detected only with favorable chopping throw and measured 
only at one epoch. For the narrow-band images, a recursive standard deviation algorithm with sigma clipping gives, 
over an empty sky area at the SE corner of our images, the following values, scaled over 1 arcsecond area: 13.1 mJy 
at 7.7 /um, 4.3 mJy at 8.7 /xm, 11.6 mJy at 9.6 yum, 6.8 mJy at 10.2 /xm, 7.5 mJy at 11.6 jim and 9.7 mJy at 12.4 fxm. 
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3. RESULTS 

The final reconstructed images are sliown in Figure 2 (10 fxm) and Figure 3 (20 /um). The size of the final 10 /xm mosaic 
is 1144 X 1072 pixels, corresponding to approximately 5.0 x 4.7 arcmin. However, as we were able to complete only the 
central columns to their originally planned extent, the rectangular area fully imaged is smaller, i.e. 1144 x 800 pixels 
= 5.0 X 3.5 arcmin. The 20 fim image has a size of 1144 x 702 pixels, corresponding to approximately 5.0 x 3.0 arcmin. 

The morphology at 10 /um is clearly dominated by the BN/KL complex and by the Ney-Allen nebula, the two 
bright sources visible at the center of the field (BN/KL to the north, Ney-Allen to the south). Extended emission 
permeates the entire region, with a general decrease of surface brightness with increasing angular distance from the 
center. This intense, diffuse emission is barely visible in the original, non-reconstructed images, since the chopping 
and nodding process filters low frequency spatial components. At 20 /im, the BN/KL complex dominates. Several 
spatially resolved knots of diffuse emission are also visible at both wavelengths. The Bright Bar, a photo-dissociation 
region viewed edge-on approximately 2' to the southwest of the Trapezium, is visible at the lower-left corner of our 
lOfiui image, in the vicinity of 9^0n A (source MAX-172). Several other filaments are visible at 10 fim in the nebular 
background. Their length range from ~ 10" to ~ 1' and their surface brightness is comparable, with peaks in the 
range ^ 0.1 — 0.2 Jy arcsec^^ at 10 /im. In the central region the peak brightness of the clumps and filaments rises to 
~ 0.7 Jy arcsec"^, in particular to the south-west and to the north-west of the Ney-AUcn nebula and to the south of 
the BN/KL, in the vicinity of the Orion-S region (see Section 4.4). 

The most remarkable new feature is represented by three, possibly four, large arcs of mid-IR emission located approx- 
imately 1 arcmin to the south of the Trapezium/Ney- Allen nebula, at R.A.(2000)=5''35™15^ DEC.(2000)=-5''24'10" 
(Figure 4a, b). They point toward the Trapezium stars (sec also Figure 12 for a sketch of their location inside the 
nebula) and appear brighter on the side facing the Trapezium. These structures are also visible at shorter wavelengths, 
in particular in the HST/WFPC2 image of O'Dell & Wong (1996), and most probably represent surfaces illuminated 
by grazing UV flux from the central stars. 

3.1. The Trapezium region 

Our 10 and 20 /um images of the Trapezium region (Figure 5a,b) show, with improved sensitivity and spatial 
resolution, the general morphology first observed at 10 /zm by McCaughrean & Gezari (1991), and studied more in 
detail by Hayward et al. (1994). The extended mid-IR emission discovered by Ney & Allen (1969) is resolved into 
a number of diffuse arc-like structures pointing toward ^^^Ori C. The brightest one, corresponding to the brightness 
peak of the Ney-Allen nebula, is centered on the BOV star 9^0ri D, undetected in our images. In the rest of this 
paper, we shall use the name Ney-Allen nebula to indicate this feature only. Another five arcs encircle the sources 
LV1-LV5 of Laques & Vidal (1979), the first photoionized disks observed in the vicinity of ^^Ori C. All arcs are 
spatially resolved in our images, and display different geometries and surface brightnesses. A most remarkable object 
lies in the immediate vicinity of O^Ovi C. This is the circular knot at the center of Figure 5a first identified as SC3 by 
Hayward et al. (1994). O^Ori C itself is also visible at 10 /im, but not at 20 ^m, as the fainter source ~ I'.'S to the 
east of SC3. All together, the arcs appear to trace the edge of a cavity centered on 9^0y\ C, filled with diffuse mid-IR 
emission, and with surface brightness increasing to the north-east. The morphology at 20 /xm is similar, but we do not 
detect the peak that Hayward (1994) found 20" to the SW of the Ney-Allen nebula at wavelengths longer than 20 /^m. 

3.2. The BN-KL complex 
3.2.1. Morphology 

In Figure 6 we show the 10 and 20 /im enlargements of the BN-KL complex from the data taken on December 18, 
2000, whereas in Figure 7 the same region is shown in the six narrow-band silicate filters. All images are absolutely 
calibrated and presented in logarithmic scale to enhance the faintest emission across the field. A few artifacts are 
present, generally related to the extraordinary brightness of the BN peak: in the broad N-band image (Figure 6a) 
there are two ghost images of BN, approximately at (4.5", —0.5") and (0.5", 6.0"), caused by inner reflections in the 
camera optics (sec also Figure 8a) . Bandwidth limitations of the readout electronics create in all narrow-band images 
an anomalous "spike," approximately 4 pixels wide, protruding to the West of BN. Some extra noise is also visible in 
the highly stretched N and Q-band images as a horizontal periodic pattern and in the vertical direction, in the rows 
and columns corresponding to the BN peak. In general, however, the images have excellent quality and are diffraction 
hmited at all wavelengths (A/D on UKIRT ranges from 0"41 at 7.7 /xm to 0"67 at 12.4 /im, and 1"08 at 19.9 /im). 
Due to the high surface brightness of the region, it is not easy to see the spectacular diffraction pattern of BN, with 
the exception of the first ring shown saturated in all images except the Q-band. Several diffraction rings appear in the 
differential images that we will present later in Section 4.3. To preserve the highest resolution, considering also that 
we are mostly dealing with the surface brightness of extended objects, we have not rebinned or convolved our images 
to our lowest resolution. In Figure 9 we show a color composite image of the BN/KL complex obtained combining our 
7.7 /im (blue), 12.4 /im (green) and 19.9 /im (red) images. 

Early mid-IR studies (Becklin & Neugebauer 1967; Rieke et al. 1973) resolved this region into 8 sources (IRCl-8, BN 
corresponding to IRcl). Another 10 compact sources were added by Gezari et al. (1998) to the IRc list, and recently 
Shuping et al. (2004) added five more. Our images have resolution and sensitivity comparable to those of Grenhill et al. 
(2004) and Shuping et al. (2004) but cover the entire field in 8 bands. Our narrow-band 12.4 /im image (Figure 7f), 
in particular, is comparable to Figure 1 of Shuping et al. (2004) and Figure 2 of Grenhill et al. (2004). We detect 
new features, especially in the unexplored outer regions of the complex. Like the large majority of compact sources 
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already cataloged, all new features appear extended and resolved into clumps, ridges, and arc-like structures down 
to our ~ 0'.'5 resolution limit. Instead of adding new entries to the already long list of IRc sources, we indicate the 
new sources using their offset (ARA", ADEC") relative to BN. For the 30 IRc sources and sub-components already 
classified, we will still use the list of Shuping et al. (2004). The position of the most relevant sources is indicated in 
Figure 8. 

BN, the main luminosity peak longward of 2 /im, remains unresolved at all wavelengths. Gaussian fits to the image 
profile at 10 /im show a ~ 3.2 pixel full-width at half-maximum (FWHM) both in right ascension and declination, 
corresponding to FWHM ~ 0'.'87. For an IR source with a spectrum increasing with wavelength, like BN, the 
diffraction limit in broad-band filters is approximately set by the longest wavelength passing through the filter, Xmax — 
13/im in our case. For an effective telescope diameter D = 3.7 m, it is X/D = 0'.'72. The pixel sampling and the 
mosaicing process, performed by integer pixels, add another ~ 1/2 pixel of image blur, bringing the FWHM right to 
the measured value. The BN peak sits on top of an extended clumpy plateau. Shuping et al. (2004) have already 
identified the remarkable SW arc, and we just add to the list a compact source immediately to the south at (—2", 2"). 

Source n of Lonsdale et al. (1982), coincident with a bipolar radio source resolved by Mentcn & Reid (1995, their 
source "L"), has been recently analyzed in detail by Grenhill et al. (2004) and Shuping et al. (2004). Using LWS on 
Keck, both team find that source n is extended, with a FWHM of approximately (0'.'75 x O'.'SO) ± O'.'IO and major axis 
oriented ~ 120°. Our pixel scale, approximately 3 times larger than that of LWS, does not sample adequately the core 
of the diffraction figure at 7.7 /im. A two dimensional Gaussian fit at this wavelength gives for source n a FWHM of 
(0'.'47 X 0'.'51), but it must be remarked that the Airy function is poorly reproduced by a Gaussian curve. If we use the 
BN peak as a reference for the PSF core and compare horizontal and vertical cuts across both sources, then source n 
remains unresolved at all wavelengths in our images. 

Other two sources, IRcl2 and IRc7 appear point-like at 7.7 /im but become extended at longer wavelengths. Sources 
IRcll and IRc6, on the other hand, show compact cores at 12.4 /xm. The point source, MAX-69 (see Section 3.4), is 
clearly visible at (11",— 6.5"), ~ 20" south-east of BN. This is a well known double star which remains unresolved in 
our images and most probably lies in the foreground of the BN/KL complex. 

Diffuse emission becomes brighter to the south-east with longer wavelengths. At 7.7 /im, sources IRc2, IRc4, and 
IRc5 appear aligned in the northeast-southwest direction and abruptly truncated to the south-cast, but at longer 
wavelengths this sharp edge, especially for IRc4, becomes less prominent and a broad plateau of emission raises around 
IRc8 and IRc22. At 19.9 /xm, extended emission fills the entire field and two dark lanes become clearly visible, one 
oriented NS along a line roughly passing through positions (8", 0") and (3", —18"), the other, almost perpendicular, 
immediately to the south of the "equatorial belt" composed of IRcl4, Ircl2, Ircll, Irc2 and Irc7. These dark lanes 
are coincident with the ammonia emission of Wilson et al. (2000), which trace the density distribution of the Orion 
hot core (Grenhill et al. 2004, see also) . Another remarkable feature is the peculiar "eastern ring" of sources centered 
at (10", -4") and traced by IRcl6, IRc 15, Irc6E, IRc21, IRc2, IRcll, IRcl2, IRcl4, Ircl3, plus other faint emission 
around (12", —2"). The area inside the ring is mostly faint clumpy emission and contains the majority of the H2O 
masers detected by Johnston et al. (1989) in OMC-1. The similar shape of the prominent IRc3 and IRc4 sources is 
also noticeable; both are arc-like and point to the IRc2 region. The tails of their arcs arc part of another remarkable 
new feature, a group of equally spaced knots tracing an arc-like "necklace" between IR5 and IRc3. The knots have all 
comparable brightness, arc elongated and radially oriented and point to the IRc2 central region. 

Finally, IRc2 is well resolved at lO/xm (Figure 6a), and shows a morphology that matches the structure of sources 
IRc2-A to -D resolved by Chelli et al. (1984) and Dougados et al. (1993) at shorter wavelengths. Its morphology has 
been discussed in detail by Grenhill et al. (2004) and Shuping et al. (2004). We will return on this object in Section 4.3. 

3.2.2. Photometry 

In Table 4, we present the photometric data at 35 positions within the BN/KL complex. The first 30 are the sources 
listed in Table 1 of Grenhill et al. (2004), the last 5 are the three brightest knots to the southwest of IRc4 (Knot 1-3), 
the compact source south of the SW Arc (BN SW) and the star MAX-69. Photometry has been obtained using a 
circular aperture of 1'.'128 diameter, corresponding to 2 square arcsec area, but the fluxes are normalized to Jy arcsec"^. 
The only exception is the BN peak, for which we give the total flux within a synthetic aperture of 5"2. All values are 
uncorrected for the background emission, and the zero point errors listed in Section 2.2.3 represent the main source 
of uncertainty. The corresponding spectra are shown in Figure 10. The silicate absorption feature is prominent at all 
positions except for the very last object MAX-69, which is probably a foreground member of the Trapezium cluster. 
A secondary maximum at the center of the absorption feature is also generally present, with the exception of the main 
BN peak and its immediate surroundings. This emission feature sits at the center of telluric ozone absorption band 
and, especially at the lowest signal levels, may be hard to detect in less than excellent weather conditions. Having 
carefully checked our absolute calibration, and verified that our photometry for the main BN peak agrees with the 
previously published spectrophotometric data (see Section 4.3), we believe this emission feature is real. 

3.3. OMC-1 South 

The OMC-IS region, located ~ 100" south of BN/KL complex, is characterized by intense outflow activity probably 
generated by deeply embedded sources (O'Dell 2001). In Figure 11, we show the position of the point sources found 
at shorter wavelengths together with the contours of the high-velocity outflow of Rodn'guez-Franco et al. (1999b). 
Part of the field is also visible in Figure 4. The Lada et al. (2000) list of candidate protostellar sources contains five 
objects in this area, TPSC-1, TPSC-2, TPSC-16, TPSC-48, and TPSC-78. Three of them, namely TPSC-1, TPSC-2, 
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and TPSC-78 are the sources with the reddest K — L color of the entire L-band survey: the first two, undetected in 
the K-band by Lada et aL (2000), have lower limits at [K-L]> 6.09 mag and [K-L]> 4.90 mag, respectively, whereas 
TPSC-78 has [K-L]= 4.85 mag. TPSC-1 is nearly coincident with our source, MAX-61 that shows significant variability 
at lO/um, with [N] = 4.40 mag in November 1998 and [N] = 3.68 mag in December 2000 (see Section 4.5.1). Object 
TPSC-2 is MAX-64, and with [N] ~ 5.3 mag it also shows a conspicuous L — N 3.8 mag color index. 

3.4. Point sources 

In Figure 12, we map the location of the 177 point sources detected at 10 and 20 /xm. In Table 3, we present their 
cross-identification, location and 10 fiva photometry. In particular, we list their sequential running number (column 
a) and their entry number in the catalogs of Hillenbrand (1997), Hillenbrand & Carpenter (2000), and Mucnch ct al. 
(2002) (column b-d), all at shorter wavelengths; coordinates (columns e,f) are derived from our final mosaic with 
typical accuracy ~ 0'.'5, accurate enough to unambiguously identify each source. In a few cases, residual pointing 
errors at the edges of the mosaic cause larger discrepancies, listed in the last column (1). The magnitudes are given 
for the two epochs (columns g,h). For a few sources, we could derive photometry only in the final mosaic (column i). 
For 22 sources, there is no counterpart identified at shorter wavelengths. Those without counterpart can therefore be 
regarded as deeply embedded objects. A few of them, however, are close to our sensitivity limits and appear only on 
the final mosaic. As explained is Section 2.2.1, we cannot rule out in these cases the possibility that they are artifacts 
of our data reconstruction algorithm. Further observations are needed to confirm them as real sources. 

4. DISCUSSION 
4.1. Diffuse emission 

With the exception of the bright BN/KL complex and the Ney- Allen nebula, the morphology of the Orion Nebula 

in the mid-IR is more similar to the radio continuum and dereddened line images obtained by O'Dell & Yusef-Zadeh 
(2000) than to the optical or near-IR images. This because the foreground "veil" of neutral material immediately in 
front of the nebula, which affects the morphology of the nebula in the visible, becomes optically thin at mid-IR (and 
radio) wavelengths, whereas the stellar emission, prominent at shorter wavelengths, becomes fainter. For this reason, 
the Dark Bay, a region of high optical depth {Ay > 3) approximately parabolic in shape, that is protruding from the 
east of the nebula up to w 30" east of the Trapezium, is essentially invisible in our images, confirming its nature as 
an enhancement of dust column density within the foreground neutral lid (O'Dell & Yusef-Zadeh 2000). 

The filamentary structures seen in emission at 10 /xm appear to be associated to the extinction-corrected color images 
of O'Dell & Yusef-Zadeh (2000). In particular, the "East Bright Bar" and "E-W Bright Bar" shown in Figure 5 of 
O'Dell & Yusef-Zadeh (2000) have 10 /zm counterparts. On the other hand, their "North Bright Bar," which is close 
to the top edge of our 10 /xm field, remains undetected, as well as the comet-like structure approximately 1' to the 
east of the Trapezium, but they are clearly visible on the L band images of both Lada et al. (2000) and McCaughrean 
et al. (1996). Most of the 10 /im filaments appear aligned along arcs roughly centered on the Trapezium. This result 
suggests that the filaments trace the edges of the main HII cavity rather than, e.g., shocks produced by collimated 
outflows. The morphological similarity between the filaments and the bars seen in the radio and in extinction corrected 
line images indicates that they are also associated with the folds in the ionization front seen edge-on. Most probably, 
the 10/im emission is dominated by PAH emission in the photodissociation region. Sloan et al. (1996) have shown 
that the 11.3/im PAH feature, produced by the aromatic C H bonds when a grain is transiently heated through the 
absorption of a UV photons, dominates the mid-IR spectrum of the Orion Bar. Similarly, the great arcs visible at 
10 /im approximately 45" south of the Trapezium (Figure 4a,b), reminiscent of the pillars seen in M16, may result 
from the advance of the ionization front between the dense and clumpy medium of the surrounding molecular cloud. 

4.2. The Trapezium region 

The complex arc-like morphology of the mid-IR emission in the vicinity of O^Ori C (Figure 5a,b) has clear optical 
counterparts in recombination lines, first detected by McCuUough et al. (1995). Using [OIII] and Ha line images 
taken with the HST, Bally et al. (1998) found that similar arcs are present aroimd many proplyds within 30" from 
O^Oti C. Bally et al. (2000) also found arcs in the outskirts of the Orion Nebula, especially in the LL Ori region (which 
remains outside of our imaged field). Hayward et al. (1994) first suggested the arcs are the bow shocks produced by 
the interaction of the fast (~ 500 km s~^) wind from 9^0v\ C with the slow (~ 10 — 50 km s~^) photoevaporated flows 
emerging from the ionization front at the proplyd surfaces. This interpretation is supported by proper motion (Bally 
et al. 2000) and spectroscopic (Hcnney et al. 2002) measurements, showing that these are stationary structures. The 
wind-wind interaction originating the arcs has been studied analytically by several authors (e.g. Canto ct al. 1996, and 
reference therein), and numerically by Garcia- Arredondo et al. (2001). There arc in fact two basic configurations: the 
interaction of two spherical winds, of interest in our case, and the interaction of a spherical wind with a plane-parallel 
flow. The latter is probably the case in the outskirts of the Orion nebula, where the wind-wind interaction is occurring 
between T-Tauri winds and the plasma flowing from the main ionization front (Bally et al. 2000). 

The analytical analysis of Canto et al. (1996), which leads to an exact solution of the wind-wind interaction within 
the "thin shell" approximation, provides an approximate solution that can be easily compared with observations. The 
ratio of the wind momenta, 

M2V2 
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where Afi,2 and vi,2 are respectively the mass loss rates and terminal velocity of the two winds, governs both the 
location of the stagnation point and the shape of the arcs. The ratio between the distance from the stagnation point 
to the source inside the arc, i?o, and the separation of the two outflow sources, D, is given by 

fio_ P'/' 
D 1 + /JV2 

The shape of the bow shock is characterized by the asymptotic angle ^oo of the tail, ^oo, given by 

9oo - tan6'oo = i- p ' '^'^^ 

If the bow shock axis is tilted with respect to the line of sight by an angle, 6i, the projection effect increases both 
the apparent ratio and the apparent asymptotic angle. For large tilt angles, the entire back side of the bow shock 
becomes visible and the characteristic arc-like geometry disappears. In these cases, the bow shock will look like a 
diffuse source in the vicinity of the central star. 

Whereas the ratio Rq/D can be easily obtained from the data, it is difficult to estimate the asymptotic angle since 
the surface brightness of the arcs decreases in the tail region. For this reason, we characterize the shape of the bow 
shock through the apparent distance, -Rgo, between the central star and the arc at 90° elongation from the direction of 
^-'^Ori C. It turns out that different combinations of /? and d.^ generate morphologies with ratios Rq/D and Rc^^/D that 
vary within a relatively narrow range. This variation is shown in Figure 13, where we have plotted the combinations 
obtained for different values of (3 and 9i. The bow shock morphology, limited to the cases with (3 =0.1, 0.04, 0.004, 
and 9i = 30°, is also shown to visualize the change of shape along the region. 

In Figure 13, we also plot the values oi Rq/D and Rqo/D measured for the Ney- Allen nebula and the five LV arcs. 
In the case of LVl, LV4, and possibly LV5, the parameters are compatible with the thin shell model. But in the 
other three cases, namely LV2, LV3, and the Ney-AUen nebula, there is a clear disagreement. In correspondence to 
the observed values of Rq/D, the values of Rgo/D should be much larger, because if the two winds carry a similar 
momentum (/? ~ 1) and the stagnation point moves closer to O^Oii C, the bow shock is expected to become more 
and more open. Instead, the three arcs with stagnation point closer to ^^Ori C show a roughly circular shape, with 
Ro — -Rgo- 

The intrinsic limitations of the algebraic solution of Canto et al. (1996) based on the thin shell approximation cannot 
explain this anomaly. Garcia- Arredondo et al. (2001) have shown that, despite the simplifying assumption of stationary 
flow and momentum conservation, the analytic treatment nicely matches their detailed hydrodynamical simulation. A 
possible explanation for the disagreement may better be provided by the detailed morphology of the arcs. The LV 3 
arc, in particular, is clearly asymmetric with respect to the direction to ^'^Ori C. It becomes narrower and brighter to 
the south in the direction of SC3, with a secondary maximum evident to the NE, at ~ 90° from the ^^Ori C direction. 
The asymmetry indicates that other dynamical effects, not accounted for by the standard models, are taking plac;e. 
A relative motion of LV3 relative to O^Oii C would break the cylindrical symmetry assumed by simple wind-wind 
models. A rather modest proper motion, of the order of a few km s~^, would be enough to add a significant tangential 
component to the speed of the photoevaporated flow, which is of the order of the sound speed in a typical ionized 
gas, Cs ~ 10 km s^^. Other dynamical effects such as the interaction with the coUimated mass outflows produced by 
the central LV stars may also locally affect the morphology and surface brightness of the arcs, leading to inaccurate 
estimates of the opening angle. 

The bright arcs of the Ney-AUen nebula are peculiar also from another point of view. Let us compare the six 
brightness profiles along the cuts joining ff^Ori C to ^^Ori D and the LVl-5 sources (Figure 14). The LVl, LV3, and 
LV4 arcs show, on the side facing ^^Ori C, brightness distributions at 10 and 20 /xm very similar to each other. The LV2 
and LV5 arcs, on the other hand, are surrounded by the diffuse luminosity, mostly originated by the Ney-AUen nebula 
and by the LV4 arc, and their 10 and 20 /xm profiles cannot be easily compared. The position of their peaks, however, 
is the same at both wavelengths, suggesting that also for them the 10 and 20 fim emission may be spatially coincident. 
The fact that the LVl-5 arcs show a similar morphology at both wavelengths supports our previous hypothesis that 
the thin shell is a good approximation and that the discrepancies between models and observations can be ascribed to 
relative motions. The case of the brightest arc (Ney-AUen nebula), however, is different. Figure 14a shows that the 
side facing 6^0ii C does not show a sharp edge at 10 or at 20 /xm. Moreover, the 20/zm emission is more extended 
in the direction of O^Oti C, indicating that a colder dust component becomes brighter at larger distances from the 
central star. The intensity profiles therefore suggest that (a) the emission is not concentrated in a thin shell, and (b) 
the UV flux from ^^Ori C is not the dominant heating source. Departures from the thin shell approximation have 
theoretical groimd. In particular, Garcia- Arredondo et al. (2001) have shown that if the stellar wind passes through 
an internal reverse shock, the wind-wind interaction takes place in a sub-sonic regime. In this case, a turbulent mixing 
layer develops at the interface between the two winds, and filaments and "fingers" are produced by the entrainment 
of the fast wind within the dense subsonic material. Instead of the sharp decay typical of a thin shell, the drop of 
surface brightness will appear in this case more gradual, and extended. Among the sources in the vicinity of ^^Ori C, 
this scenario applies naturally, and exclusively, to the Ney-AUen nebula, as it is associated to the BOV star O^Ori D. 
This hot star originates a substantial stellar wind, and the formation of a internal reverse shock where the gas fiows 
enter a subsonic regime is a natural outcome in the early adiabatic phase of expansion (Castor et al. 1975; Weaver 
et al. 1977). 
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The main peculiarity of the Ney-AUen nebula, however, is that it contains in an extremely harsh environment a 
considerable amount of dust, estimated between ~ 2 x 10^'' g if = 200 K and ~ 1 x 10^® g if = 350 K (Hayward 
1994). Proplyds in Orion are supported by the photoevaporation of circumstellar disks around low mass pre-main 
sequence stars (Churchwell et al. 1987; Meaburn 1988; Henney & Arthur 1998(@; Johnstone et al. 1998; Henney & 
O'Dell 1999(@). In analogy, it is natural to think to a massive circumstellar disk around ^^Ori D as the most natural 
dust reservoir for the Ney- Allen nebula. The presence of a massive circumstellar disk is compatible with the standard 
view that high-mass stars reach the main-sequence while they are still in their protostellar accretion phase (Stabler 
et al. 2000). In this scenario, it is important to notice that O^Oii D is the only star of the Trapezium without a binary 
or multiple companion, down to a resolution limit X/D = 57 mas and 76 mas in the H and K bands, respectively 
(Preibisch et al. 1999). Conversely, no other Trapezium stars seems to be associated with excess mid-IR emission. 
This result may therefore indicate that for OB stars the presence of low mass companions is not compatible with 
substantial circumstellar disk, and vice versa. 

Concerning the SC3 source, i.e. the bright extended blob in the immediate vicinity of O^Ori C, we have discussed 
in a previous paper (Robberto et al. 2002) its apparent morphological anomaly, as it appears undisturbed by the 
powerful 9^ On C wind despite having the shortest projected distance from ^^Ori C. On the basis of its high mid-IR 
brightness, we have proposed that it may represent a disk on the far side of the Nebula, oriented face-on with respect 
to both ^^^Ori C and the Earth. High spatial resolution mid-IR observations, possibly with interferometric techniques, 
should easily resolve the radial distribution of its bright mid-IR emission and allow us to clarify the nature of this 
most enigmatic object. 

4.3. The BN/KL complex 
4.3.1. Main physical parameters 

The iibiquity and strength of the silicate absorption feature in the 10 /im region indicates that foreground absorption 
plays a dominant role in shaping the observed morphology of the BN/KL complex. As a baseline, one can therefore 
assume a two component model (Aitken et al. 1981), where a background blackbody source is shielded by a foreground 
absorbing medium too cold to contribute significantly to the mid-IR emission. This assumption allows us to readily 
obtain a map of the foreground opacity from the line to continuum ratio at 9.8 /xm. In a first order approximation, 
the continuum can be derived from the interpolation at 9.8 /im of the fluxes measured at 7.7 and 12.4 /im, whereas 
for the depth of the silicate feature we interpolate between the fluxes measured at 8.7 and 10.2 /im to avoid the 
secondary maximum seen in the 9.6 /im filter. The natural logarithm of the absorption line strength is the optical 
depth and is presented in Figure 15. The optical depth in the main central region increases from west to east across 
the entire field, peaking in the immediate vicinity of the dark lane. It reaches the maximum value at the IRc2 position, 
Tg.s = 3.37 (unless otherwise specified, all values are per pixel, i.e. averages over an area ~ 0'.'06, corresponding to 
^9.8 — 3.66. This is ~ 20% lower than the value reported by Gezari et al. (1998). The drop of extinction coinciding 
with the dark lane is due to the absence of flux sufficiently high to allow for a meaningful estimate of the extinction 
and must therefore be regarded as a selection effect. 

The map of the peak optical depth allows us to correct all images for the foreground extinction. There are several 
possible choices for the spectral profile of the silicate feature. We adopt here the "Orion Trapezium" emissivity 
(Forrest et al. 1975), as updated by the spectroscopic observations of Banner et al. (1995). We average the Banner 
et al. emissivity over the (cold) filter transmission profiles provided by the manufacturer. We also include the non- 
silicate extinction, that we parameterize in terms of the visual extinction, Ay, using a fit to the extinction data of 
Rieke & Lebofsky (1985), A\ = 0.381 • Ay ■ A~^'^^^. Following Rieke & Lebofsky, we weight the two contributions by 
assigning to the silicate extinction at the 9.8 /im peak a 90% fraction of the total extinction at that wavelength. 

The morphology of the BN/KL complex after extinction correction for all 6 silicate filters is shown in Figure 16. 
As expected, the changes with respect to the unreddened images (Figure 7) are stronger where the silicate absorption 
feature is more prominent. In particular, at 9.8 and 10.2 /xm, the IRc2 object recovers its rank of second brightest 
source at all wavelengths after BN peak. In Figure 17, we zoom in on the IRc2 region, comparing for all 6 silicate 
filters the morphology before and after extinction correction. The most striking result is the appearance of a point 
source just to the south of the IRc2 peak. It was visible in the raw 7.7 /im images, but now becomes evident also in 
the dereddened 9.6 and 10.2 /im images. Its position is coincident with knot IRc2-A of Dougados et al. (1993) and lies 
~ 0'.'5 north of radio source I of Menten & Reid (1995), which remains undetected. The fact that in the dereddened 
images the mid-IR peak is visible at 7.7 /xm, invisible at 8.7 /um, visible again in the dereddened 9.6 and 10.2 fj,m 
images may be due to a combination of imperfect extinction correction, intrinsic silicate emission and loss of resolution 
at the longest wavelengths. 

From the dereddened images, we derive the distribution of the 10 /im color temperature of the underlying emission 
by fitting, pixel by pixel, a blackbody curve to all 6 silicate filters. The result is shown in Figure 18. The temperature 
across the BN/KL complex typically ranges between 180 K and 300 K, with maxima at the BN peak with Tcoi — 1000 K, 
and at the IRc2 position with Tcoi — 420 K. The BN temperature is certainly overestimated. Since BN is unresolved 
and our images are diffraction limited, the encircled energy decreases with wavelength. The dilution of flux over a 
larger number of pixels causes an apparent increase of the color temperature at the peak. Concerning IRc2, the warm 
dust appears mostly elongated in the northeastern-southwestern direction, i.e. in a direction parallel to the dark lane, 
and peaks ~ 0'.'5 to the southeast of source I. This effect is shown in Figure 18, where we overplot the contours of the 
dereddened 7.7 /im image at 1.5 and 3.0 Jy arcsec"^ flux levels. 
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The bolomctric luminosity can finally be estimated from the L = crT'^^i • iiTR^ relation, where R is the radius 
corresponding to 1 resolution element (0.26" ~ 120 AU at the distance of the Orion Nebula). We do not show its 
spatial distribution, which is very similar to Figure 18 with an even more prominent role of IRc2 and BN. 

In Table 5, we list the values of the three parameters we have derived — optical depth at 9.6 /xm, color temperature 
in the 7.7-12.4 /xm band and total luminosity — for the sources with photometry presented in Table 4. The values of 
Tg.e and Tcoi are averages over a circular aperture of 1 arcsecond square area, whereas the luminosity, L, is integrated 
over the same area. Table 5 also reports the temperatures and luminosities derived from the dereddened 10 and 20 /im 
broad-band filters having assumed at 20 /um the interstellar extinction law of Mathis (1990). The color temperature 
between these two wavelengths (125-140 K) is typically 80 K lower than that derived from the silicate filters (200- 
230 K). The discrepancies up to a factor of 10 in the estimated luminosities reflect these difi^erences. 

The elementary technique we have used so far, based on a crude estimate of the extinction, still provides maps of 
the spatial variations that are quite robust, since the data have been taken all at the same time and calibration errors 
are systematic. In the next section we shall discuss a more refined model to investigate further the characteristics of 
the mid-IR emission of the BN/KL complex. 

4.3.2. Silicate absorption feature 

The two component model has been successfully used in the past to model the 10 /im spectrum of the main BN 
peak (Aitken et al. 1981, 1985). Figure 19 compares our broad-band photometry of BN calculated over a synthetic 
aperture of 5'.'2 diameter with the spectrophotometric data of Aitken et al. (1985) and the synthetic photometry of 
Gezari et al. (1998). The agreement between different data sets is good, considering that the broad-band values have 
not been color corrected and the differences in aperture diameter between different authors. The solid and dotted 
lines represent the locus of a blackbody curve at Tcoi before and after reddening by an optical depth rg.g to fit our 
photometry. We use the Trapezium emissivity of Banner et al. (1995) with the additional A~^-^^ continuum discussed 
above. The best fit to our data has T^oi = 373 K and rg.s = 1-37, with a filling factor FF = 5.41 x 10~^ to normalize 
the absolute flux levels. From the filling factor, we derive a photospheric radius Rbn = 27.2 AU at 450 pc, and the 
resulting luminosity is Lbn = 47r_R^^ • aT^^i = 600 Lq. For comparison, Aitken et al. (1981) derive Tcoi — 336 K and 
Tg.s = 1-54, whereas the total luminosity of BN is Lbn — (1 — 2) x W^Lq (Scoville et al. 1983)). 

It is known that, besides BN peak, the two component fit does not account for the detailed shape of the 10 /xm spec- 
trum of the other compact sources of the BN/KL complex (Aitken et al. 1985). The 8-13 /xm absorption feature 
appears to change across the complex, both in depth and in shape (Aitken et al. 1981). These variations reflect differ- 
ent physical conditions in the absorbing medium. Downes et al. (1981) and Werner et al. (1983) first indicated that 
most of the IRc sources discovered in the BN/KL complex are not self-luminous, but rather extended clumps of diffuse 
emission. Genzel & Stutzki (1989) indicated that the BN/KL complex is dominated by the radiative transport of IR 
radiation in a highly non-uniform medium. A further complication may arise from the fact that at shorter wavelengths 
scattering becomes important. Werner et al. (1983) were the first to report that at 4.8 /im the BN/KL complex is 
predominantly a refiection nebula, but since the albedo of astronomical silicates drops by two orders of magnitudes 
between 4.5 and 8 /xm (Draine 1985), one can still assume that scattering is negligible at mid-IR wavelengths. In this 
case, the basic two component model can be modified to account for the following effects: 

1. The background emission is normally optically thin. Only for a few sources, like BN and IRc2, the background 
emission may be represented by as a stellar photosphere multiplied by a proper filling factor. 

2. The presence of foreground absorbing medium (as in the two component model). 

3. A third component, optically thin, needed to account for silicate emission at 9.6 /xm. 
In a slab geometry, these three factors are described by the equation: 



F^^FF*B^{Tb) 1 - exp-^'^'^) exp"^"^'') -hSj.(T,) 1 - exp 



(4) 



Here FF < 1 and are the filling factor and temperature of a background medium with optical depth, t;,. For a 
photosphere (e.g. BN-peak), it is [1 — exp"'^''] ~ 1; Ta is the optical depth of the cold foreground absorbing material, 
that we assume is mostly located in the vicinity of the background sources (i.e. part of the Orion molecular cloud); 
Te and Te are the temperature and optical depth of the foreground gas responsible for the silicate emission. We locate 
this third component in the foreground of the main absorbing layer, assuming that it physically corresponds to warm 
dust dust heated in the foreground cavity, i.e. contained in the optically visible Orion Nebula. 

Equation (4) contains six free parameters, and the choice of different emissivity laws allows for some extra freedom. 
For this exploration, we adopt the silicate emissivity law of Jager et al. (1994), which provides excellent agreement 
with the spectra of massive young stellar objects and the Ney Allen nebula in what concerns the position of the 10 /im 
peak and has a relatively narrow width that may provide the best agreement with the most recent measures of the 
Trapezium (see below). Figure 19 shows the best fit to the flux (averaged over 4 pixels) at the IRc2 peak. To allow 
the numerical routine to converge, we had to reduce the degree of complexity of our system. Thus, the dotted line 
represents the fit to the 7.7, 8.7, 10.2, 11.6 and 12.4 bands, i.e. all the silicate filters except the 9.6 /im one, provided by 
the first term at the right-hand side of Equation (4), having assumed FF = 1. We obtain = 710 K, Tb — 0.0019, and 
To = 21.6. The dashed line represents the contribution from the second term at the right-hand side, having assumed 
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ad-hoc Ts = 300 K and Tg — 1.5 x 10^"^ to approximately match the photometry in the central part of the silicate 
feature. The thick solid lines represent the sum of the two contributions. Figure 19 shows that the 20 fim point cannot 
be matched by a single blackbody temperature, due to the strength of the secondary silicate absorption peak at 18 /zm. 
At 10 /xm, the emissivity law of Jager et al. is too broad to match the 8.7 and 10.2 points, but the addition of 
the third component in emission provides a reasonable fit, including the 9.6 /itm point. A detailed look shows that the 
emissivity law of Jager et al. (1994) is still too broad to reproduce at the same time the 9.6 /um peak and the adjacent 
8.7 and 10.2 fiia points. The finite filter band-pass, together with the presence of a range of temperatures instead 
of the assumption of a single 300 K blackbody, exacerbates the problem as both effects tend to reduce the observed 
sharpness of the feature. To illustrate the effect of the bandpass, we have convolved the sum of the two contributions 
with a AA = 1 fj,m rectangular bandpass. The result is presented as a thin solid line in Figure 19. Similar results 
are found for other sources, and in Figure 19b we show as an example the case of IRc7, where we find Tf, = 232 K, 
n = 0.04, Ta = 14.2, and = 0.0035. 

The value rg.g — 22 we have found for lRc2 corresponds to ~ 400 and is compatible with the extreme value 
of the optical depth r > 300 at 8.0 /xm at the near ammonia core estimated by Grcnhill et al. (2004). A general 
emissivity law with a narrower peak would provide higher values of column density. On the other hand, there is no 
reason to assume a-priori that the emissivity is the same for the three components, or in general across the region. 
There is evidence that, for yomig stellar objects and star forming regions, the 10 fiia silicate feature depends on the 
environment, due to the interplay between dust column density, cosmic ray irradiation and grain coagulation (Bowey 
et al. 2003). In particular, the shape of the canonical "Trapezium emissivity" of Forrest et al. (1975); Roche & Aitken 
(1984); Whittet et al. (1988), generally assumed as a reference for molecular clouds and used to constrain the optical 
constants of the popular Draine & Lee (1984) model, shows a dependence on the spatial scale. Unpublished arcsecond 
resolution spectra of the Trapezium obtained by T. Hayward (Bowey 1998) show that the peak of the 10 /im feature 
is significantly narrower (0.5-0.6 iim) than that reported by Forrest et al. (1975) and Banner et al. (1995), these last 
authors with a 9'.'4 beam (beam size is unquoted in the Forrest et al. paper). Different broadening may be due to grain 
size effects or to changes in the grain structure or optical properties (see Jager et al. 1994, for a discussion). Bowey 
(1998) suggests that measures made with large apertin-c, sampling a wide range of dust temperatures, may also cause a 
broadening of the 10 /im silicate feature, at least in the case of Trapezium. If the optically thin emission is due to dust 
grains exposed to the UV radiation, the grain structure may be affected in the direction of a predominantly crystalline 
composition. Amorphous silicates may be converted to crystalline minerals when annealed (heated) in vacuum at high 
temperature (Day & Donn 1978). The Trapezium grains are also exposed to strong UV radiation, but since they may 
have been only recently released by the disk surrounding ^^-D, they may still represent a relatively standard mixture of 
amorphous and crystalline silicates. On the other hand the dust grains inside the Orion nebula, and those distributed 
in front of the BN region in particular, could have been exposed for a longer time to the nebular UV radiation. It is 
intriguing to notice that the reddening law in Orion has Ry = Ay/E{B — ]/) ~ 5.5 vs. the standard value 7?v' = 3.1 
of the interstellar medium. This anomaly could be related to the peculiar evolution of the dust grains in the harsh 
nebular environment. 

4.3.3. The role of IRc2 and source n 

IRc2 has been considered for a long time as the mid-IR source more closely associated to a deeply embedded center 
of activity of the region, e. g. the origin of the spectacular outflow phenomena emerging from the BN-KL complex. 
Evidence in this direction comes from the position of IRc2 at the center of (1) an expanding cluster of maser sources 
with velocity centroid nearly coincident to the LSR velocity of the hot core (Genzel et al. 1981); (2) a well defined 
polarization pattern (Aitken et al. 1997; Chrisostomou et al. 2000), and (3) a rich outflow phenomenology observed at 
ncar-IR wavelengths (Johnes & Walker 1985). Source I of Menten & Reid (1995), a compact HII region approximately 
1" south of IRc2, has been more recently indicated as the actual center of activity of the BN/KL complex (see 
Tan 2004, for a recent review). Whether IRc2 is just a complex of knots heated by source I or, as Grenhill et al. 
(2004) speculate, a group of individual stars associated with source I in a young multiple system, remains unclear. 
An intermediate case could also be considered, with source I and the IRc2-A peak representing deeply embedded 
individual sources, and the other IRc2 knots being diffuse sources heated mostly by IRc2-A peak. 

Source n, whose radio emission has been resolved into a bipolar nebula by Menten & Reid (1995), has also recently 
gained a prominent role. It is only 3"away from IRc2, lies exactly at the center of expansion of the H2O masers of 
Genzel et al. (1981), and within the narrow 2'.'0 x 0'.'5 strip where Gaume et al. (1998) locate their "shell" masers. 
Both Shuping et al. (2004) and Grenhill et al. (2004) have been able to resolve this source at 10 nm with Keck. It 
appears slightly elongated in a direction perpendicular to the axis of the bipolar radio source. Shuping et al. (2004) 
suggest that source n lies at the center of a bipolar cavity excavated by the same outflow responsible for the reach 
maser phenomenology. The presence of a compact HII region clearly suggest that this is a young early type star. 
Grenhill et al. (2004) estimate a luminosity L ^ 2000 Lq (mid-B type), but also notice that the recent detection of 
hard-X ray emission is unusual for early- type stars. 

Gczari et al. (1998) have used mid-IR photometric data similar to ours to downplay the role of IRc2 vs. source n 
as the dominant energy source of the region. In particular, Gezari et al. (1998) estimate for IRc2 a total luminosity 
L ~ 1O^L0, approximately 2 orders of magnitude lower than the previous estimates (e.g. Downes et al. 1981; Wynn- 
Williams et al. 1984). One may question the reliability of a technique that provides values different by orders of 
magnitudes in front of a rather modest gain in spatial resolution. In fact, these estimates must always be taken with 
great caution. First, we have seen in the previous paragraphs that different analysis techniques may lead to very 
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difFcrcnt values for the optical depth at 9.8 /im (rg.S = 3.37 vs. Ta = 21.6), even if both use the same data set (in our 
particular case, both neglect the 9.6 /im point). It is also known that the commonly used relation, L = AirR^aT^e, 
where R and T are the source radius and temperature and e represents the dust emissivity, may lead to a significant 
underestimate of the bolometric luminosity if one assumes the IR color temperature as the blackbody temperature of 
the source (Panagia 1975). A more subtle effect may also be present, since the same relation also assumes that the 
size of the underlying photosphere does not change with wavelength. In the case of IRc2, however, if the absorbing 
medium is clumpy or in general inhomogeneous, an unresolved source may have a filling factor that increases with 
wavelength as the foreground absorbing material becomes more and more transparent. The underlying source will then 
appear redder, and the estimated luminosity will be lower due to the strong dependence on the color temperature. In 
particular, a value Ay = 60"* corresponds to A9.7 ~ 3.4"* at 9.7 /xm and to A20 = 1-3'" at 20 /um, using Bessel & Brett 
(1988) IR reddening law. Thus, IRc2 can be nearly optically thick at 10 /zm and nearly optically thin at 20 /xm. In 
these conditions, the mid-IR morphology is determined by a complex combination of thermal emission and variations 
of foreground extinction. Further multiwavelength observations will need to be coupled with 3-d radiative transfer 
models to disentangle each contribution. 

4.4. OMC-1 South 

The Orion-S region is characterized by the rich phenomenology typically associated with massive star formation: 
warm gas (Ziurys et al. 1981), large dust column densities (Keene et al. 1982), H2O masers (Gaume et al. 1998), 
molecular outflows (Schmid-Burgk et al. 1990; Rodriguez-Franco et al. 1999a), and several optical jets and Herbig 
Haro objects (Bally et al. 2000). With 9 mid-IR sources in less than 1' x 1' field (Figure 11), this is a rather dense 
area and the scale of the outflows, up to a projected distance of ~ 100" = 0.21 pc in the case of IIII203/204 often 
makes it difficult to unambigously associate each optical outflows to a well defined driving source. We have already 
mentioned the correspondence between the reddest 3.5 fim sources, TPSC-1 and TPSC-2 of Lada et al. (2000), with 
MAX-61 and MAX-64, respectively. A system of outflows that includes HH 529 and HH 269, and possibly HH 528, 
appears to depart from an "optical outflow source" (named OOS by O'Dell & Doi 2003a) coincident, within the error 
ellipse, with both sources. Based on their proper motions, O'Dell & Doi (2003a,b) estimate for the OOS coordinates 
RA(2000.0)=5''35'"14!56, DEC(2000)=-5°23'54", within ~ 3" from the position of MAX-61 and MAX-64'l MAX-61, 
which is radio source 144-351 of Zapata et al. (2004), is the most viable candidate, since our images (Figure 4) show 
two elongated filaments immediately to the east and west of MAX-61 with the same orientation of the HH 529/HH 269 
outflows. These arc most probably the jets emerging from the source, with the mid-IR emission due to entrained dust 
on their outside, along with gas. Doi et al. (2004) combine proper motions with new radial velocities of the outflows 
to determine that the optical jets originate from a source lying only about 0.02 pc within the PDR. Their relations 
between proper motion, radial velocities and the mid-IR emission would provide a self- consistent picture of this activity 
center, except that both systems of HH objects appear blue-shifted and breaking through the ionization front, a rather 
unusual circumstance that may indicate a more complex scenario. A curvature of the outflow toward our line of sight 
could result from the interaction of the jet with some stellar wind, originated e.g. from the Trapezium stars or from 
the nebular ionization front. This is the hypothesis of Smith et al. (2004). It is clear that if the jet direction can 
be deflected by wind interaction, the identification of candidate outflow sources on the basis of proper motion and 
jet/counterjet alignment becomes more uncertain. In fact. Smith et al. (2004) cast doubts about the association of 
HH 269 with HH 529 and rather suggest that HH 269 is the counter lobe of the HH 203/204 pair, with significant 
deflection with respect to both the plane of the sky and our vantage point. 

Concerning the other prominent optical outflows, both Zapata et al. (2004) and Smith et al. (2004) argue that 
HH 202 is associated to HH 528 rather than HH 203/204, as suggested by Rosado et al. (2001). Zapata et al. (2004) 
suggest that this system originates from their radio source 143-353. At this position, we detect an extended knot of 
10 emission shown saturated in Figure 4 and not included in our list of 177 point sources. Smith et al. (2004) 
suggest a different scenario, with the HH 202 jet associated to their source nr. 1 (MAX-45), but they also note that 
neither MAX-61 nor MAX-64 cannot be ruled out. 

In what concerns the molecular outflows, the A = 1.3 mm bolometer maps of Mezger et al. (1990) show a peak 
of dust continuum brightness, source FIR4, located at RA(2000)=5''35'"13!3, DEC(2000)=-5°24'12'.'5. Low velocity 
molecular outflows depart from the vicinity of FIR4 (Ziurys et al. 1990; Schmid-Burgk et al. 1990). A second, very 
young (dynamical age of ^ lO'' years), very fast 110 km s^^) and very compact (< 0.16 pc) molecular outflow has 
also been discovcircid by (Rodrfguez-Franco et al. 1999a, b). The geometry of the high velocity gas suggests that its 
powering source is located ~ 20" north of FIR4. We detect no IR sources within a 12" radius from FIR4, in particular 
there is no mid-IR source at the position of the 1.3 cm source 134-411 of Zapata et al. (2004). At the other position 
20" north, there are three near-IR sources of Hillenbrand & Carpenter (2000), i.e. HC 190, HC 198 and HC 199, but 
none of them have been detected at longer wavelengths either by us or by Lada et al. (2000). On the other hand, 
sources MAX-42 and MAX-43, corresponding to TPSC-16 and TPSC-46 of Lada et al. (2000) and to the radio sources 
136-359 and 136-356 of Zapata et al. (2004), are displaced ~ 5" — 10" to the south-east, still in a favorable position 
with respect to the symmetry point of the CO outflow. Both show a remarkable ii" — L ~ 3 color and appear as 
excellent candidates for the exciting source. The fainter of the two, MAX-42, is better aligned with the symmetry axis 
of the outflow. The best ground-based near-IR images clearly show MAX-43 is diffuse. Source MAX-47, corresponding 

^ Zapata et al. (2004) base their comparison with the present paper on a preliminary version of the manuscript. In this final version 
their sources MAX-58 and MAX-60 are now Usted as MAX-61 and MAX-64 
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to radio source 139-357 of Zapata et al. (2004), also lies along the axis of the outflow. We have detected this source 
only at 20 /im, whereas Smith et al. (2004) have only a 10 fiia detection. Since there is no near-IR counterpart at this 
location, this appears to be a deeply embedded source and therefore also represents a good candidate for the source 
driving the molecular outflow. 

The last protostellar candidate of Lada et al. (2000), TPSC-78, corresponding to source MAX-46, is located near a 
characteristic feature of the Orion Nebula, the dark arc limiting to the north the Orion-S region, « 45" to the south- 
west of the Trapezium. Since the dark arc also appears in the 20 cm continuum image, it is not due to foregromid 
extinction but rather traces a curved part of the ionization front seen edge-on, similar to the linear Orion Bar (O'Dell 
& Yusef-Zadeh 2000). It is remarkable to find to the north of dark arc (i.e. on the "wrong side") one of the brightest 
{[N] = 2.74"*) and reddest (L — ~ 4.3, N — Q ~ 2.1) sources of our survey. This object can still be embedded 
within the high density region creating the dark arc, if one assumes that the arc traces the closest part of a neutral 
cloud, whose surface is not exactly edge-on but rather tilted with respect to our vantage point. As we have mentioned 
above. Smith et al. (2004) consider this bright object as the most likely driving source for the HH 202/HH 528 outflow 
system. 

4.5. Point sources 
4.5.1. Mid-IR variability 

For several 10 fim sources, the comparison of the photometric results obtained in November 1998 versus December 
2000 suggests the presence of mid-IR variability. In Figure 21 we plot against each other the data for the two epochs, 
listed in Table 3. Whereas the average magnitude difference between the two runs is less that 0.02 mag, the average 
scatter is ~ 0.27"*, much larger than the average measure error A[N] ~ 0.11"'. This discrepancy is too large to be 
explained by measurement errors unaccounted for. In fact, for this type of comparison our photometric errors are 
probably overestimated, as they include, added in quadrature, the uncertainties on the absolute zero point. These are 
systematic shifts common to all magnitudes taken on the same night, that we conservatively estimated of the order 
of 0.05-0.07 mag at 10 /im. Since almost all the data presented in Figure 21 were taken in two nights, November 
27, 1998 and January 18, 2000, and since the systematic differences appear to be less than 0.02 mag, our zero points 
uncertainties may actually be too generous, hardening the evidence for variability. 

Stellar variability should be common among these sources. Typical members of the Trapezium clusters are spotted, 
rotating T Tauri stars, known to be variable at visible wavelengths. According to Herbst et al. (2002), essentially all 
stars in Orion Nebula bright enough to be monitored from the ground are variable, with periods having a bimodal 
distribution with peaks at 2 and 8 days for M> O.25M0, and unimodal distribution peaked near 2 days for M< O.25M0. 
Future comparative studies of the optical and mid-IR variability, based on simultaneous multiwavelength photometric 
monitoring, may allow to understand whether the variability is due to disk activity, e.g. fluctuations of the mass 
accretion rates, or to changes in the foreground extinction due to inhomogeneities of the circumstellar environment or 
of the circumestellar disk itself. 

4.5.2. Mid-IR excess 

To investigate the nature of the mid-IR emission of our sources, we compare in Figure 22 the [K-L] vs. [N-Q] color 
indexes of our Trapezium stars versus a sample of classical T Tauri stars in the Taurus- Auriga association. For the 
Trapezium stars, we complement our photometry with the K and L-band photometry of Muench et al. (2002), using 
two different symbols to indicate the 10 nm data obtained in the 1998 and 2000 observing runs. 

The Taurus- Auriga sources are those listed in Table II of Strom et al. (1989). Their K and L-band photometry is the 
average of the values available for each fllter in the recent literature, mined using the VizieR Service^, and in general 
turns out very close to the photometry of Strom et al. (1989). The 10 and 20 /xm photometry is derived from the IRAS 
photometry at 12 and 25 /xm in the following way. The IRAS PSC2 fluxes are quoted at the effective wavelengths 
assuming an input energy distribution constant per log v, i.e. i>F,^ = constant. This is also the shape of the typical 
"flat" Spectral Energy Distribution (SED) of circumstellar disks aroimd T Tauri stars. We therefore assume the IRAS 
fluxes to be color corrected and adjust their values only to account for the difference in effective wavelength between 
the filters. From the i/fy=constant condition, we derive m(10) = rn(12) -|- 0.516 and m(20) = rn(25) -|- 0.694 using the 
zero magnitude fluxes reported in Section 2 of this paper and in the IRAS PSC2 manual. We have visually inspected 
each SED plotting, together with all the IR data available in the literature, our interpolated 10 and 20 /im points. 
Whereas our points are always in excellent agreement with the overall shape of the SED (admittedly constrained by 
the IRAS points), most of the 10 and 20 /im ground based data are old and suspiciously away from the best fit line. 

Figure 22 clearly shows that there is a systematic difference between the Orion and the Taurus-Auriga samples, as 
the Orion stars have a [N-Q] excess between 3 and 5 magnitudes, whereas for Taurus- Auriga the [N-Q] index ranges 
between 1 and 3 magnitudes. Also, the [K-L] excesses appear on average higher for Orion, although for both samples 
the majority of stars have [K-L] in the range 0.7-1.5 mag. We show in Figure 22 also five stars without L-band 
photometry, but well detected at shorter wavelengths, to which we have assigned [K-L]=0. Also for these sources the 
[N-Q] index appears higher than that of the stars in Taurus- Auriga. 

Extinction cannot explain these effects, either in terms of intensity or reddening. In Figure 22 we have indicated 
with an arrow the reddening corresponding to Ay = 10 magnitudes of extinction. Hillenbrand (1997) suggests that 

® VizieR Service is a joint effort of CDS and ESA-ESRIN, available on the internet at http://vizier.u-strasbourg.fr/viz-bin/VizieR 
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the mean extinction for the opticaUy visible population of the Orion Nebula cluster is Ay ^ 2.5 magnitudes, with few 
optical stars having Ay > 5 magnitudes. For the IR sample, the mean extinction is approximately 3 times higher, but 
only a small fraction of sources has Ay > 15 magnitudes (Hillenbrand & Hartmann 1998). 

Calvet et al. (1992) have calculated the IR colors of accreting disks around classic T Tauri stars assuming a variety of 
stellar and disk parameters. Each one of their models is represented in Figure 22 as a solid line, as different orientation 
angles are taken into account. The Calvet et al. (1992) predictions fall in general within the region of the observed 
T Tauri star colors. In particular, the models with the highest [N-Q] > 2 have in common a higher value of 7, the 
index of the power law describing the increase of the disk height with the distance from the central star. Large disk 
flaring angles could therefore explain the extreme [N-Q] colors of our sources. 

In a previous paper (Robberto et al. 2002), we have explored the continuum SED of circumstellar disks embedded 
within a HII region like Orion. We predicted that the extra energy input from the ionizing stars and from the 
nebular environment raises the disk temperature and causes an increase of the disk flaring angle. An inflated disk 
will also intercept a larger fraction of the flux emitted by its own central star. In the vicinity of ^^^Ori C, the neutral 
flow produced by disk photoevaporation originates a neutral envelope bordered by a ionization front, the "proplyd" 
phenomenon (O'Dell et al. 1993; O'Dell & Wen 1994). Together with the disk, the dusty envelope also contributes 
to the mid-IR emission. Our models show that there is in general a strong enhancement of the mid-IR emission due 
both to the higher disk temperature and the presence of the photocvaporated envelope. The UV heated dust produces 
a characteristic far-IR peak of emission at 30-60 /im. In Figure 22 we show as a dotted line the locus of the colors 
resulting from our model, assuming the limiting case of a disk facing both the ionizing star (for simplicity O^Oii C 
is considered the only ionizing source) and the Earth at various distances from ^^Ori C. The [N-Q] colors predicted 
by this model nicely match those of the Orion sample, supporting the hypothesis that the UV radiation is ultimately 
responsible for the strong enhancement of 20 /im emission. 

In Robberto et al. (2002), we have shown that the IR SEDs of proplyds depend on a number of factors, both physical, 
like the distance from the ionizing star, the size of the envelope, the inner disk radius, the temperature of the central 
star and the dust composition, and geometrical, like the tilt angle of the disk with respect to the ionizing star and to 
the Earth. The projection factors in particular have a dramatic impact on the observed SED, since disks surfaces at the 
same distance from 9^0n C can be exposed to very different UV fluxes and appear from the Earth under all possible 
angles, showing either the front or the back side. The 30-60 /xm peak stands as the main observational characteristics 
of these systems. Depending on the dust temperature and grain size, IR color indexes that do not include the 20 /xm 
band, like [K-N] or [L-N], may remain unaffected by a mid-IR peak too cold to produce strong excess at 10 /xm. For the 
main group of sources detected only at 10 /xm, it is necessary to analyze in detail each individual source to constrain 
the disk properties. This analysis will be presented in a second paper, and will include the result of our deep 10 /um 
survey of selected fields in Orion. 

5. CONCLUSION 

Using a new restoration method developed for this project, we have obtained the first wide-field and diffraction 
limited 10 fim and 20 fim images of the core of the Orion Nebula. The images confirm that the morphology of the 
Orion Nebula at mid-IR wavelengths is dominated by the BN/KL complex and by the Ney- Allen nebula. At lower 
signal levels, the Orion Nebula shows a filamentary structure that traces the walls of the HII region. A remarkable 
group of arc-like structures « 1' is detected to the South of the Trapezium. We analyze the structure of the Ney-AUen 
nebula in the framework of analytic wind-wind interaction models, finding evidence for a complex kinematical status 
at the center of the Cluster. The B0.5V star 6'-'^0ri-D, associated with the brightest mid-IR structure of the Trapezium 
region, is probably surrounded by a photoevaporated circumstellar disk. The fact that, among the main four Trapezium 
stars, this star is the only one without a binary companion suggests that stellar multiplicity and massive circumstellar 
disks are mutually exclusive in the Trapezium. We detect in the OMC-1 South region three point sources that can be 
associated to the mass loss activity detected at millimeter wavelengths. Finally, we list the position and photometry 
of 177 point sources, the large majority of which detected for the first time in the mid-IR. Twenty two of them lack 
counterparts at shorter wavelengths, and are, therefore, candidates for deeply embedded protostars. The comparison 
of photometric data obtained at two different epochs reveals that source variability at 10 /xm is present up to ~ 1™ 
level on a time-scale ~ 2 yr. We compare the mid-IR color-color diagram for our sources in Orion and a similar dataset 
in the Taurus association, finding indications for large mid-IR excess, especially at 20 /zm, and discuss this finding in 
the framework of the Robberto et al. (2002) model. 

We have also presented narrow band images of the BN/KL complex taken in the standard silicate filter set. The 
BN/KL complex is resolved into multiple and extended features. In particular, we resolve a point source in cor- 
respondence of the IRc2-A knot. We derive the extinction and temperature maps on the basis of the simple "two 
component" model for the IR emission, and discuss the shape of the silicate absorption feature and possible scenarios 
for dust evolution in the HII region. 

The authors are indebted to the MPIA and UKIRT staff for their support with the MAX operations. The contribution 
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Fig. 1. — Example of the construction of the 10 fim mosaic. Clockwise from top-left, individual chopped and nodded frames are 
progressively assembled into a final strip. The restoration method works on each indivial column (or row). The column, approximately 35' 
wide, contains the BN/KL complex and the OMG-IS region. North is at the top and west is on the right. 
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Flc;. 2. — Final 10 /im mosaic of the Orion nebula. The two saturated sources at the center of the field are 1) to the southeast, the 
Ney-AUen nebula, in correspondence of the Trapezium stars, and 2) to the northwest, the BN/KL complex. 
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Fig. 3. — Same as Figure 1, at 20 fxra. The only prominent source at this wavelength is the BN/KL complex. 
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Fig. 4. — Arcs of mid-IR emission approximately 45" south of the Trapezium (see Figure 11 for their location). Left: original image 
before image reconstruction. Right: same field after image reconstruction and filtering. Scale is 45" x 45", with standard orientation. The 
sources in the north-western corner fall within the OMC-1 South region (see Figure 11). 




Fig. 5. — The Trapezium region at 10 fj,m (left) and at 20 /im (right). O^Ovi C and SC3 are the two circular sources at the center of the 
10 /im image, O^Ori C being the fainter one to the left. Field size is 1' x 1', north is up and east is to the left. 
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Fig. 6. — The BN/KL complex at 10 /im (left) and 20 /im (right). The scale is logarithmic and ranges between 0.1 and 10 Jy/arcsec^ at 
10 //m, and between 3 and 300 Jy/arcsec^ at 20 /um. 
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Fig. 7. — The BN/KL complex at 7.7 and 8.7 (top row), 9.6 and 10.2 fiva {middle row), and 11.6 and 12.4 fim {bottom row). All 
images are shown in logarithmic scale and units are ,Iy arcsec^'^. To facilitate the comparison, all contours are plotted at the same levels of 
0.1, 0.25, 0.50, 0.75, 1, 2.5, 5, 7.5, and 10 Jy arcsec"^. Equatorial coordinates are relative to the peak of the BN source (see Section 2.2.2). 
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Fig. 8. — Contour plot of the BN/KL complex at 10 fim (left) and 20 (right), with labels indicating the main features discussed in 
the text. Contours levels are 0.22, 0.30, 0.37, 0.45, 0.60, 0.80, 1.0, 1.2, 1.35, 1.5, 1.87, 2.25, 3.0, 4.5, 6.0, 10, 15, 30, 60, 100, 150 and 200 
Jy/arcsec2 at 10 fim, and 7, 10, 12, 14, 16, 18, 20, 25, 30, 35, 40, 50, 60, 80, 100, 120, 140,1 60, 180, 200 Jy/arcsec^ at 20 fim. 




Fig. 9. — Color composite of the BN/KL region obtained combining the 7.7 ^m (blue), 12.4 fim (green) and 19.9 ^m (red) images. 
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Fig. 10. — Mid-IR spectral energy distribution for the 35 sources of Table 4. Errors refer to the typical zero point error of our absolute 
calibration. The wide horizontal line indicates the broad-band N filter. 
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14^ 13^ 
RA(2000.0)=5^35'^ 

Fig. 11. — Map of the sources located in the Orion-S region, together with the intensity contours of the CO outflow mapped by Rodriguez- 
Franco et al. (1999b). The naming conventions are: 1) MAX sources: this paper; 2) TPSC sources: Lada et al. (2000); 3) HC sources: 
Hillenbrand & Carpenter (2000); 134-411 (et simiJia) and DOS sources: Zapata et al. (2004); FIR 4: Rodn'guez-Franco et al. (1999a). 




Fig. 12. — Location of our mid-IR point sources. Thin open circles: 10 fim sources with counterpart at shorter wavelengths. Thick open 
circles: 10 fira sources detected at 10 fim and not at shorter wavelengths. Filled circles: sources detected only at 20 fim. Small circles: a 
subset of the near IR sources that remained undetected at 10 fim, presented here for reference and ease of comparison with near-IR images. 
Thin solid lines trace the main filamentary structures and the arcs discussed in Section 4. B^Ori C is source 111; S^Ori A is source 178. 
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Fig. 13. — Locus of the morphological parameters Rq/D and Rgo/D describing the Trapezium arcs. Open circles refer from left to right 
to /3 = 0.001, 0.002, 0.004, 0.01, 0.02, and 0.1, and increase from the bottom to the top from 6i = 0° up to the maximum apparent tilt angle 
in steps of 15°. The location of the arcs in the vicinity of O^Ovi C is also shown, with the observational uncertainties. 
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Fig. 14. — Brightness profiles along the cuts joining 6^0ri C to 6^0ri D and the LVl-5 sources. Solid line: 10 fira emission; dashed line: 
20 /Ltm emission, scaled for comparison. 
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Fig. 15. — Map of the optical depth of the BN/KL complext at 9.6 //m. 
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Fig. 16. — Same as Figure 7, after reddening correction. 
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Fig. 17. — Comparison of the morphology of the IRc2/IRc7 region before and after reddening correction for the six silicate filters. 
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Fig. 18. — Map of the 10 fim color temperature 
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Fig. 19. — Silicate spectrum of BN. Diamonds: this paper; Plus signs: spectrophotometry of Aitken et al. (1985); asterisks: broad-band 
data of Gezari et al. (1998). The solid line refer to a blackbody curve at T^-^i = 373 K. The dotted line joins the blackbody points 
corresponding to our effective wavelengths after extinction correction with optical depth rg. 8=1-37. 
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Fig. 20. — Fit to the mid-IR photometry of IRc2 (left) and IRc7 (right) based on Equation 4. Dotted Hne: first term on the right side of 
Equation 4, relative to the bagground source and foreground extinction. Dashed line: scncond term on the right side of Equation 4, relative 
to the optically thin emission of 300 K dust. Thick solid line: sum of the two contributions. Thin solid line: sum of the two contribution 
convolved with a 1 /im rectangular bandpass to reproduce the effect of finite bandwith. Error bars represent the zero point errors quoted 
in Section 2.2.3 
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Fig. 21. — Comparison of the MAX [N]-band photometry at the two epochs. 



28 



5 - 



at 



1 



' I ' ' ' M 

X 



X + 



+ 
+ 




I I . I I . I 



I . . i . I 



-0,5 0,0 0,5 1,0 1,5 2.0 2.5 3,0 

Fig. 22. — Mid-IR color-color plot for the Orion and Taurus-Auriga sources. The K- and L-band data for Orion are from Muench et al. 
(2002). The cross (x) and plus {+) signs refer to the MAX 10 fim data taken in 1998 and 2000, respectively. Filled circles are for the 
Taurus dataset. Solid lines represent the colora of the Calvet et al. (1992) flaring disk models. The diamonds and the relative dotted line 
represent the model of Robberto et al. (2004) for a photoevaporated disk in a HII region oriented face-on with respect to both the ionizing 
star and the observer. The points refer to different distances from the ionizing star, from top to bottom: 30", 60", 90", 120" and 150", 
corresponding to 0.065, 0.13, 0.20, 0.26 and 0.33 pc, respectively, at the assumed distance d = 450 pc of the Orion Nebula. 



TABLE 1 
Main observational parameters 



Night 


Filter 


I.T. 


Chop. Preq. 


I.T./frame 


Chop. Throw 






(ms) 


Hz 


s 


arcsec 


27 November 1998 


N 


20.2 


4.2 


20.07 


9.80 (N-S) 


28 November 1998 


N 


20.2 


4.2 


40.14 


18.55 (N-S) 


18 December 2000 


N 


9.9 


3.9 


25.36 


24.64 (E-W) 


19 December 2000 


N 


9.9 


3.9 


25.36 


24.64 (E-W) 


19 December 2000 


Q 


9.9 


3.9 


12.68 


24.64 (E-W) 


18 December 2000 


silicates 


20.2 


4.12 


12.944 


24.64 (E-W) 



TABLE 2 

Calibrators for narrow-band imaging 



Source 


7.7 /im 


8.7 //m 


9.8 fj,m 


10.2 ^im 


11.2 //m 


12.4 ^lm 


HR 1040 


2.55 


2.55 


2.55 


2.55 


2.55 


2.55 


HR 1457 


-2.97 


-2.94 


-3.00 


-3.04 


-3.07 


-3.07 


HR 1713 


0.06 


0.06 


0.06 


0.06 


0.06 


0.06 


HR 4069 


-0.84 


-0.83 


-0.91 


-0.96 


-1.00 


-1.00 


HR 2943 


-0.77 


-0.77 


-0.76 


-0.76 


-0.76 


-0.76 


HR 5340 


-3.11 


-3.10 


-3.14 


-3.17 


-3.18 


-3.18 


F^=o(Jy) 


63.56 


50.86 


41.64 


36.67 


28.94 


25.28 



TABLE 3 

MiD-IR Photometry of Sources in the Orion Nebula 



MAX 


H97 


HCOO 


MLLAOO 


R.A. 


Dccl. 


[N] 


[N] 


[N] 


[Q] 


Note 










(J2000.0) 


(J2000.0) 


(Nov. 1998) 


(Dec. 2000) 


(mosaic) 




1 


265 




234 


5''35'"05!52 


-5°24'09'.'9 




6.55±0.14 






MAX source is +3.53"E,-0.28"N 


2 


266 




223 


5'"35'"05!55 


-5°24'14'.'5 




7.50±0.39 






MAX source is +3.44"E,-0.19'N 


3 


274 




314 


5''35'"05?74 


-5°23'44y5 




6.72±0.16 






MAX source is +2.10"E,0.19"S 


4 


278 




161 


5''35'"05!74 


-5°24'18''4 




5.90±0.11 






MAX source is +1.55"E,0.32"S 


5 


291 


99 


161 


5''35'"06!45 


-5°24'39yo 






6.85±0.15 






6 


341 


305 


409 


5'*35'"10!23 


-5°23'21'.'0 




4.70±0.09 








7 

8 




320 


425 


5''35'"iar23 
5''35'"10!28 


-5°21'43'.'5 
-5°23'15'.'8 




6.48±0.18 
7.46±0.14 








9 


349 


146 


220 


5'' 35"* 10! 56 


-5°24'17'.'2 




6.78±0.10 








10 


347 


554 


733 


5''35'"lffi64 


-5°21'56'.'8 




6.16±0.13 








11 


392 


385 


507 


5'' 35"" 10? 64 


-5°22'55'.'8 




6.14±0.09 








12 


352 


224 


316 


5''35'"10?75 


-5°23'44'.'l 




7.27±0.13 








13 


355 


416 


542 


5''35'"10?90 


-5°22'46'.'l 




6.37±0.10 








14 


356 


93 


148 


5/135771 ;^ OP 94 


-5°24'48'.'5 




6.00±0.12 








15 




415 


541 


5''35'"11^04 


-5°22'46''3 




6.50±0.11 








16 








5'\35'"11!11 


-5°21'53"0 




6.38±0.17 








17 




434 


575 


5'*35'"11!20 


-5°22'37'!'6 




5.06±0.09 








18 




105 


171 


5'"35'"11!32 


-5°24'38y4 




6.34±0.13 








19 




203 


289 


5'"35'"11!49 


-5°23'51'.'5 




6.95±0.10 








20 


366 


140 


207 


5''35'"11!65 


-5°24'21''9 




7.99±0.19 








21 




611 


809 


5''35'"11!75 


-5°21'31"8 






7.23±0.3 






22 




672 


880 


5'*35'"11!88 


-5°20'59"9 






5.70±0.2 






23 


3058 


664 


872 


5''35'"11?92 


-5°21'02'.'8 






6.85±0.3 






24 






710 


5''35'"11!95 


-5°22'01''5 






7.17±0.4 






25 


9009 


395 


516 


5'"35'"11?98 


-5°22'53{'6 




6.03±0.10 








26 


373 




944 


5''35'"12?01 


-5°20'33'.'8 






6.10it0.3 






27 


378a,b 


215 


302A,B 


5''35'"12!27 


-5°23'47'.'8 




6.00±0.10 








28 


377a,b 


698 


914A,B 


5'' 35"" 12! 30 


-5°20'45'.'6 






5.24±0.3 






29 




761 


793 


5/135771 ]^2!41 


-5°21'36'.'7 




8.08±0.25 








30 


9018 


368 


482 


5'' 35"" 12! 58 


-5°23'02'.'0 




7.43±0.15 








31 


385 


228 


318 


5''35'"12!60 


-5°23'43'.'9 




5.49±0.10 




1.39±0.20 




32 


391 


699 


916 


5/135771 ]^2!85 


-5°20'44'.'8 






8.97±0.3 






33 




606,608 


804A,B 


5/135777 ;^2!85 


-5°21'33'.'6 


6.33±0.13 










34 


9029 


505 


673 


5''35'"13!07 


-5° 22' 16'.' 1 


7.36±0.31 


7.38±0.29 








35 


399a,b 


487 


651A,B 


5''35'"13!19 


-5°22'21'.'8 


7.27±0.19 


7.04±0.18 








36 








5''35'"13!19 


-5°22'04y6 


7.14±0.12 


7.35±0.19 








37 


401 


80 


136 


5''35'"13!20 


-5°24'55^'5 




3.79±0.11 








38 


400 


682 


896 


5''35'"13!24 


-5°20'54"7 


6.34±0.20 










39 








5'"35'"13!31 


-5°21'53{'6 


7.55±0.26 










40 




273 


369 


5'*35'"13!40 


-5°23'29"5 


6.22±0.09 


6.21±0.11 








41 


409 


240 


327 


5 35 13? 44 


— 5 23 40. 2 


D.44±0.U9 


D.22±0.10 








42 




178 


263 


5'"35'"13!55 


-5°24'00^'2 


6.74±0.09 


6.87±0.12 








43 




192 


276 


5''35'"13!58 


-5°23'55"0 


6.40±0.09 


6.24±0.11 




2.9±0.3 




44 




602 


797 


5'>35'"13!71 


-5°21'35'.'7 


6.43±0.09 


6.39±0.13 








45 


424 


541 


715 


5''35'"13!79 


-5°22'03'.'5 




6.28±0.50 








46 


431 


242 


328 


5''35'"13!79 


-5°23'40'.'l 


2.69±0.08 


2.78±0.09 




0.68±0.19 




47 








5'»35'"13!8 


-5°23'57'.'0 








1.55±0.18 




48 


423 


703 


703 


5''35'"13!81 


-5°22'07'.'8 


4.90±0.08 


4.70±0.09 









TABLE 3 — Continued 



MAX 


H97 


HCOO 


MLLAOO 


R.A. 
(J2000.0) 


Decl. 
(J2000.0) 


[N] 

(Nov. 1998) 


[N] 

(Dec. 2000) 


49 




548 


726 


5''35'"13!81 


-5°22'00''6 


1.98±0.08 


1.90±0.09 


50 








5''35'"13!9 


-5°22'12'.'0 






51 


9053 


629 


828 


5''35'"13!96 


-5°21'22'.'8 


5.99±0.08 




52 




314 


413 


5'*35'"13!97 


-5° 23' 20'' 1 






53 




134 


199 


5''35'"14!01 


-5°24'25('4 




7.47±0.21 


54 








5''35'"14?1 


-5°23'00'.'0 






55 




704 


707 


5'"35'"14!12 


-5°22'05'.'4 


6.78±0.50 




56 








5''35'"14!18 


-5°23'14'.'8 




5.99±0.21 


57 


437 


135 


203 


5/135171 ]^4S29 


-5°24'24'.'7 


5.57±0.09 


5.66±0.10 


58 




361 


476 


5'' 35"" 14? 30 


-5°23'04'.'3 


6.04±0.09 


6.51±0.11 


59 


9061 


345 


452 


5''35'"14?31 


-5°23'08'.'4 


5.81±0.08 


5.51±0.10 


60 


9062 


399 


518 


5/s 35m 14333 


-5°22'53'.'7 


6.16±0.09 


6.18±0.14 


61 






292 


5''35'"14?38 


-5°23'50'.'6 


4.40±0.08 


3.68±0.09 


62 


441 


258 


350 


5'"35'"11M0 


-5°23'33''8 


5.83±0.08 


6.05±0.11 


63 


9069 


369 


485 


5'\35'"11?58 


-5°23'02''6 




6.97±0.19 


64 




189,193 


273A,B 


5''35'"14!54 


-5°23'55'.'7 


5.44±0.08 


5.27±0.09 


65 




411 


534 


5'^35'"14!63 


-5°22'50''2 




6.95±0.19 


66 


9073,9075 


714 


467A,B 


5'^35'"14!84 


-5°23'05"0 


6.58±0.19 




67 


445a,b 




922A,B 


5''35'"14!69 


-5°20'43''7 






68 


9071 


657 


863 


5'"35'"14!71 


-5°21'05{'9 


6.32±0.10 




69 


454 


431 


568 


5''35'"14?91 


-5°22'39''6 


4.79±0.08 


4.75±0.10 


70 




437 


580 


5''35'"15!19 


-5°22'37'.'3 


6.59±0.10 


6.53±0.26 


71 


463 


398 


514 


5''35'"15!21 


-5°22'54''7 


4.44±0.08 


4.28±0.10 


72 








5'^35'"15!34 


-5°22'50"4 


7.21±0.17 




73 


470 


504 


671 


5''35"»15?34 


-5°22'16'.'3 


4.85±0.08 


4.60±0.10 


74 


472a,b-9096 


475,476 


630A,B,C 


5/135171]^ 5335 


-5°22'25'.'7 


7.00±0.11 


7.07±0.32 


75 


473 


646 


848 


5'' 35"" 15? 40 


-5°21'13'.'6 


5.45±0.08 




76 


475 


223 


311 


5''35'"15?43 


-5°23'45'.'0 


6.81±0.09 


7.26±0.16 


77 








5''35'"15?49 


-5°23'52'.'3 


7.31±0.17 




78 








5'' 35"" 15? 50 


-5°23'48'.'2 




6.21±0.11 


79 


5178 


251 


340 


5/135777 1^5852 


-5°23'37'.'5 


6.51±0.09 


5.97±0.12 


80 








5/735777^585 


-5°23'10'.'3 






81 


480 


172 


249 


5''35'"15!64 


-5°24'03''9 


6.26±0.09 


6.49±0.11 


82 


479 


389 


503 


5'"35'"15!65 


-5°22'56'.'7 


3.93±0.08 


3.87±0.09 


83 


488w 


307 


405 


5'\35'"15!72 


-5°23'22'.'7 


5.75±0.09 


5.47±0.11 


84 


492 


420 


545 


5'\35'"15!78 


-5°22'46'.'4 


8.33±0.82 




85 


1865 


336 


440 


5''35'"15!80 


-5°23'14'.'6 


5.43±0.25 




86 


489 


287 


381 


5''35'"15!80 


-5°23'26''7 


5.44±0.08 


5.28±0.11 


87 


488a 


306 


404 


5''35'"15?86 


-5°23'22'.'5 


4.48±0.08 


4.46±0.11 


88 




291 


385 


5'"35'"15!86 


-5°23'25"0 


5.65±0.18 




89 


490 




925 


5''35'"15!88 


-5°20'39''7 






90 


9118 


145 


213 


5'' 35™ 15? 90 


— 5°24'18"0 




6.91±0.12 


91 


496 


490 


653 


5''35'"15?93 


-5°22'21'.'8 


6.26±0.09 




92 


499 


213 


296 


5''35'"15?95 


-5°23'49'.'9 


3.18±0.08 


3.40±0.09 


93 








5''35'"15?97 


-5° 22' 36'.' 1 


6.86±0.22 




94 


503 


202 


288 


5/73577716801 


-5°23'52'.'8 


5.16±0.08 


5.08±0.10 


95 








5''35'"16?05 


-5°22'33'.'4 


6.68±0.10 




96 


9128 


393 


511 


5''35'"16?07 


-5°22'54'.'4 


5.50±0.08 


5.07±0.10 


97 


1863b 


350 


460B 


5''35'"16?08 


-5°23'07'.'2 


3.63±0.09 


3.66±0.09 



[N] 
(mosaic) 



[Q] 



Note 



6.73±0.3 



-0.64±0.16 
0.98±0.22 



1.65±0.21 



0.42±0.18 



1.91±0.22 



MAX source is -2.80E,+0.29" S 

MAX source is -2.25"E,-0.52"N 
MAX source is -3.35"E,0.24"N 



6.60±0.3 



8.83±0.3 



1.39±0.20 
1.59±0.21 



1.41±0.24d 



1.41±0.24d 



1.42±0.19 



1.44±0.19 



OW-158-326^ 



TABLE 3 — Continued 



MAX 


H97 


HCOO 


MLLAOO 


R.A. 
(J2000.0) 


Decl. 
(J2000.0) 


[N] 

(Nov. 1998) 


[N] 

(Dec. 2000) 


[Nl 
(mosaic) 


L^J 


Note 


98 




296 


389 


5''35'"16?09 


-5°23'24'.'2 


3.93±0.08 


3.98±0.10 




1.18±0.19 




99 




504a,b 


926 


5''35'"16!13 


-5°20'38'.'4 






8.13±0.3 






100 








5''35'"16!15 


-5°22'26'.'2 


6.86±0.19 








to be confirmed 


101 


1863a 


354 


460A 


5''35'"16!15 


-5°23'07'.'l 


5.76±0.37 


6.52±0.10 








102 


511a,b 


522 


689A,B 


5''35'"16?28 


-5°22'10{'5 


5.63±0.08 


5.65±0.09 








103 


9135 


652 


857 


5'*35'"16';28 


-5°21'08'.'9 


6.40±0.10 










104 


513 


412 


535 


5'"35'"16!31 


-5°22'49'.'0 


6.62±0.33 


6.97±0.19 








105 


512 


322 


422 


5''35'"16!32 


-5°23'16'.'7 


4.24±0.08 


4.35±0.09 




1.06±0.26 




106 






401B 


5'' 35"" 16? 34 


-5°23'22'.'6 


3.06±0.08 


3.07±0.09 




0.27±0.17 


SC3'' 


107 








5''35'"16?4 


-5°22'56'.'0 








1.24±0.20 




108 


515 


171 


246 


5''35'"16?41 


-5°24'04'.'2 


4.94±0.08 


5.15±0.09 








109 




112 


173 


5'>35'"16?41 


-5°24'36'.'9 






8.10±0.3 






110 




514 


682 


5'' 35'" 16? 42 


-5°22'13'.'0 


5.52±0.08 


4.88±0.09 








111 


1891 


309 


401A 


5''35'"16?46 


-5°23'22'.'7 


4.23±0.08 


4.62±0.09 








112 








5''35'"16!51 


-5°22'48'.'0 




6.79±0.65 








113 








5'"35'"16!51 


-5°23'30y2 


6.07±0.31 








to be confirmed 


114 




325 


426 


5'"35'"16!52 


-5°23'15y5 


5.45±0.25 










115 








5'"35'"16?7 


-5°22'57yo 








1.12±0.21 




116 


524 


323 


424 


5'"35'"16!76 


-5°23'16y5 


3.86±0.08 


4.26±0.09 




0.46±0.23 




117 








5'"35'"16!82 


-5°23'58('2 




5.97±0.74 








118 




289 


383A,B 


5'"35'"16!84 


-5°23'26^'0 


2.15±0.08 


2.29±0.10 




-0.64±0.18 




119 


532 


414 


537 


5''35'"16!98 


-5°22'49^'l 




6.95±0.21 








120 








5''35'"16!98 


-5°22'31"0 




6.97±0.42 








121 


534 


252 


341 


5'"35'"16!99 


-5°23'37yi 


4.90±0.08 


5.40±0.09 




0.89±0.21 




122 


533 


375 


489 


5''35'"16?99 


-5°23'01'.'3 


6.32±0.09 


5.94±0.12 








123 


537 


244 


330 


5''35'"17?03 


-5°23'40'.'l 




6.65±0.19 








124 


538 


259 


348 


5''35'"17?06 


-5°23'34'.'0 


4.45±0.08 


4.24±0.10 








125 








5''35'"17?11 


-5°23'26'.'6 




5.09±0.12 




0.78±0.25 




126 








5''35'"17?11 


-5°23'29'.'3 




5.02±0.11 






to be confirmed 


127 


1889 




423A,B 


5''35'"17?26 


-5°23'17'.'0 




5.97±0.46 








128 


548 


440 


584 


5''35'"17?37 


-5°22'35'.'6 


5.52±0.08 


5.43±0.10 








129 


549 


180 


261 


5''35'"17?39 


-5°24'00'.'l 


6.03±0.08 


e.OOitO.lO 








130 


9180 


313 


411 


5'"35'"17!46 


-5°23'2iyi 


5.41±0.09 


4.89±0.34 








131 


554 


295 


387 


5'^35'"17!55 


-5°23'25'.'0 


5.94±0.09 


5.92±0.11 




1.7±0.3 




132 


553a,b 


388 


499.2A,B': 


5'^35'"17?57 


-5°22'57'.'10 


5.67±0.09 


5.72±0.10 








133 


552a,b 


586 


773A,B 


5'^35'"17!57 


-5°21'45'.'9 


5.30±0.08 


5.15±0.11 








134 




563 


743 


5'"35'"17!62 


-5°21'54y3 


6.33±0.09 


6.69±0.17 








135 


558 


241 


325A 


5''35'"17!68 


-5°23'4iyi 


6.44±0.10 










136 


9195 


234 


320 


5''35'"ir;75 


-5°23'42{'7 


6.88±0.11 










137 


561 


681 


892 


5'*35'"17!82 


-5°20'54{'2 






7.77it0.3 






138 


9201 


367 


480 


5'"35'"17!90 


-5°23'02('0 


5.91±0.32 










139 




542 


714 


5''35'"17!90 


-5°22'03y4 


6.16±0.09 


6.07±0.11 








140 


567 


425 


549 


5''35'"17?92 


-5°22'44'.'6 


6.13±0.09 


6.27±0.11 








141 


9196 


230 


317 


5''35'"17?93 


-5°23'44'.'5 


7.46±0.34 










142 




679 


890 


5''35'"17?96 


-5°20'55'.'8 






6.34±0.3 






143 


9210 


174 


248 


5'»35'"18!04 


-5°24'03'.'0 


4.86±0.08 


4.76±0.10 








144 


575 


177 


256 


5''35'"18?08 


-5°24'01'.'2 


6.43±0.12 


6.53±0.10 








145 


MS- 107 


372 


486 


5''35'"18ni 


-5°23'02'.'8 


6.71±0.23 










146 


581 


253 


344 


5''35'"18?21 


-5°23'36'.'0 


4.00±0.08 


4.44±0.09 









TABLE 3 — Continued 



MAX 


H97 


HCOO 


MLLAOO 


R.A. 
(J2000.0) 


Decl. 
(J2000.0) 


[N] 

(Nov. 1998) 


[N] 

(Dec. 2000) 


[Nl 
(mosaic) 


[Q] Note 


147 




533 


705 


5''35'"18!24 


-5°22'05'.'7 


7.18±0.15 


7.05±0.17 






148 








5''35'"18!31 


-5°22'22'.'6 


7.38±0.16 








149 


9220 


348 


455 


5''35'"18!36 


-5°23'06'.'7 


6.71±0.17 








150 


582 


331 


431 


5''35'"18!37 


-5°23'15'.'7 


6.93±0.20 








151 








5''35'"18?65 


-5°23'43^'l 




7.85±0.17 






152 


596 


337 


441 


5'*35'"18?65 


-5°23'13^'9 


3.28±0.08 


3.45±0.09 






153 


598a,b 


713 


499.1A,B"= 


5'"35'"18!70 


-5°22'56'.'9 


5.16±0.08 


5.37±0.11 






154 


607 


214 


298 


5''35'"19!04 


-5°23'49'.'6 




3.91±0.09 






155 




351 




5''35'"19?6 


-5°23'07'.'O 








2.8 ±0.25 


156 


622 


133 


193 


5''35'"19?67 


-5°24'26'.'5 




5.71±0.10 






157 


624 


492 


648 


5''35'"19!84 


-5°22'21'.'9 




6.32±0.10 






158 


3075 


346 


450 


5'>35'"20?31 


-5°23'06'.'9 






7.57±0.3 




159 


648a,b 


277,279 


366A,B 


5'' 35"" 20? 45 


-5°23'29'.'9 




6.21±0.10 






160 


9266 


308 


398 


5'"35'"20M8 


-5°23'22"9 




7.60±0.18 






161 


652 


141 


208 


5'\35'"20?56 


-5°24'20''9 




7.37±0.14 






162 


658 


97 


154 


5'"35'"20!62 


-5°24'46'.'l 




5.78±0.17 






163 


660 


590 


777 


5'"35'"20!74 


-5°21'45'.'2 




6.23±0.13 






164 


661 


560 


737 


5'"35'"20!80 


-5°21'56'.'l 




6.64±0.14 






165 


665 


574 


753 


5'"35'"20!95 


-5°21'51'.'5 




6.02±0.12 






166 


669 


217 


300 


5'"35'"21!02 


-5°23'49'.'2 




7.51±0.18 






167 








5'"35'"21!33 


-5°23'42y7 




6.86±0.13 






168 


687a,b 


585,255 


769A,B 


5'"35'"21?67 


-5°21'48('l 




5.86±0.13 






169 


694 


249 


332 


5'"35'"21!78 


-5°23'39"5 




7.12±0.12 






170 


698 


204 


282 


5'"35'"21!79 


-5°23'53'.'8 




7.57±0.21 






171 


9294 


358 


465 


5''35"»21?84 


-5°23'06'.'9 






8.01±0.3 




172 








5''35'"21?94 


-5°22'28'.'9 




7.28±0.42 






173 


5179 


121 


182 


5/135171 22? 10 


-5°24'32'.'9 




6.48±0.10 






174 


707 


139 


201 


5''35'"22?20 


-5°24'25'.'0 




6.71±0.11 






175 


710 


155 


224 


5'' 35"" 2232 


-5°24'14'.'5 




6.66±0.11 






176 


721 


233 


319 


5'' 35"" 2254 


-5°23'43'.'8 






7.90±0.3 




177 


723 


339 


443 


5''35'"22?80 


-5°23'13'.'8 




5.82±0.09 






178 


1993 


82 


129A,B 


5'' 35"" 22? 90 


-5°24'57'.'5 




4.97±0.18 






179 


744 


262 


347 


5'"35'"23!78 


-5°23'34y6 






6.28±0.3 




180 


9317 


460 


601 


5'"35'"24!41 


-5°22'31'.'9 






6.41±0.3 




181 


756 


110 


164 


5'"35'"24?44 


-5°24'41'.'0 




5.50±0.21 






182 


762 


119 


176 


5''35'"24!70 


-5°24'36'.'0 




6.31±0.31 






183 


766a,b 


387 


494A,B 


5''35'"25!07 


-5°22'59'.'4 






6.01±0.3 





b From O'Dell & Wen (1994) 
^ From Hayward ct al. (1994) 

The identification number 499 is used tlie list of Muench et al. (2002) for two different stars (both are binaries). 

Source nr. 84 and 87 are unresolved in the Q-band 



TABLE 4 

Peak Strface Briciitxess of the BN/KL complex 



Source JName 


Ortset rl.A. 


Unset JJec. 






Narrow-band 






Broad-band 


7.7 


8.7 


9.8 Atm 


10.2 fim 


11 O , m\ 


12.4 fiiii 


N-band 


Q-band 


BN 


0.0 


0.0 


122.30 


106.11 


62.59 


57.21 


126.80 


143.19 


119.58 


171.34 


BN SW arc 


-1.4 


-0.6 


9.82 


11.35 


7.22 


10.31 


16.29 


20.78 


8.65 


59.76 


IR n 


3.3 


-9.8 


2.11 


1.06 


1.62 


1.22 


3.23 


7.14 


2.12 


31.53 


IRc2 A+B 


5.7 


-7.0 


15.21 


1.07 


1.38 


0.89 


5.62 


20.24 


3.52 


34.07 


IRc2 C 


4.6 


-7.8 


7.98 


0.73 


1.50 


0.95 


4.67 


14.94 


2.77 


36.67 


IRc2 D 


4.8 


-7.0 


10.37 


1.58 


1.59 


1.14 


6.67 


19.64 


3.85 


41.70 


IRc2 E 


4.5 


-6.3 


4.55 


1.08 


1.02 


0.86 


3.96 


10.17 


1.86 


28.45 


IRc3 N 


-1.6 


-7.2 


0.62 


0.44 


2.32 


1.64 


4.02 


7.75 


1.94 


66.69 


IRc3 S 


-1.4 


-7.9 


0.69 


0.43 


2.71 


1.75 


4.32 


8.81 


2.23 


74.34 


IRc4 


1.1 


-11.1 


1.63 


0.83 


3.58 


2.72 


10.13 


20.15 


4.76 


101.36 


IRc5 


0.5 


-14.5 


0.56 


0.30 


2.08 


1.30 


3.45 


7.56 


1.73 


59.99 


IRc6 


0.4 


-4.5 


0.91 


0.42 


2.13 


1.40 


3.92 


7.74 


2.11 


54.24 


IRc6 N 


-0.2 


-3.2 


0.89 


0.61 


2.30 


1.64 


3.58 


5.97 


1.99 


52.92 


IRc6 E 


1.8 


-4.8 


0.75 


0.24 


1.43 


0.96 


2.81 


6.13 


1.52 


38.26 


IRc7 


2.8 


-7.8 


4.29 


1.21 


2.84 


1.89 


8.86 


23.88 


4.93 


73.93 


IRc8 


8.4 


-12.2 


0.13 


0.07 


1.02 


0.63 


1.07 


2.04 


0.58 


28.24 


IRcll 


8.5 


-7.4 


1.38 


0.32 


0.89 


0.64 


1.90 


5.28 


1.19 


20.76 


IRcl2 


9.9 


-7.9 


1.58 


0.22 


0.72 


0.51 


1.45 


3.64 


0.85 


16.94 


IRcl3 


14.1 


-4.8 


0.14 


0.06 


0.56 


0.35 


0.55 


1.14 


0.29 


14.37 


IRcl4 


13.8 


-6.7 


0.18 


0.09 


0.78 


0.49 


0.96 


1.87 


0.52 


19.42 


IRcl5 


3.1 


-2.3 


0.45 


0.23 


1.30 


0.87 


1.91 


3.46 


1.15 


33.98 


IRcl6 N 


4.9 


1.3 


0.71 


0.24 


1.17 


0.81 


1.83 


3.09 


0.99 


30.66 


IRcl6 S 


5.0 


-0.3 


0.51 


0.18 


1.17 


0.77 


1.77 


3.26 


1.37 


31.12 


IRcl6 W 


4.0 


0.2 


0.47 


0.23 


1.20 


0.83 


1.61 


2.73 


1.58 


30.73 


IRcl7 


3.0 


2.3 


0.42 


0.26 


1.17 


0.77 


1.35 


2.08 


0.69 


28.38 


IRcl9 


5.2 


4.7 


0.74 


0.24 


0.68 


0.44 


0.97 


1.37 


0.56 


16.06 


IRc20 


0.7 


-6.9 


0.74 


0.31 


1.36 


0.90 


2.22 


5.17 


1.10 


37.28 


IRc21 


4.4 


-5.0 


0.84 


0.19 


0.64 


0.46 


1.29 


3.21 


0.67 


16.82 


IRc22 


6.4 


-11.3 


0.23 


0.08 


0.72 


0.46 


0.76 


1.83 


0.50 


17.30 


IRc23 


6.3 


2.2 


0.38 


0.12 


0.56 


0.35 


0.84 


1.59 


0.45 


15.81 


Knot 1 


-4.2 


-10.4 


0.17 


0.07 


0.96 


0.53 


0.82 


1.68 


0.42 


25.13 


Knot 2 


-3.9 


-12.7 


0.21 


0.08 


0.82 


0.48 


0.66 


1.53 


0.38 


23.44 


Knot 3 


-2.6 


-14.6 


0.25 


0.10 


0.95 


0.54 


0.81 


1.80 


0.45 


25.70 


BN SW 


-2.0 


-2.0 


1.21 


0.88 


1.90 


1.54 


3.12 


5.05 


1.90 


36.66 


MAX-69 


11.4 


-16.6 


0.13 


0.15 


0.65 


0.42 


0.40 


0.48 


0.38 


11.17 



TABLE 5 

Physical parameterh of the IRr sources in the BN/KL complex 



7.7-12.4 10-20 urn 



Smirfp IVamp 


'^9.6 


^ col 




1 ; 

col 








(K) 




(K) 




BN 


0.51 


490 


9.21 xl05 


296 


7.98x10* 


BN SW arc 


0.33 


288 


8.48x10* 


176 


8.39x10^ 


IR n 


1.30 


227 


2.41x10* 


149 


4.37x10^ 


IRc2 A+B 


2.91 


308 


1.14xlOS 


164 


6.27x103 


IRc2 C 


2.53 


258 


4.49x10* 


153 


4.70x103 


IRc2 D 


2.36 


219 


2.11x10* 


160 


5.68xl03 


IRc2 E 


1.94 


208 


1.63x10* 


148 


4.06x103 


IRc3 N 


1.11 


204 


1.46x10* 


127 


2.27xl03 


IRc3 S 


1.18 


209 


1.61x10* 


128 


2.35xl03 


IRc4 


1.53 


184 


9.55x103 


138 


3.12x103 


IRc5 


1.32 


206 


1.53x10* 


126 


2.15x103 


IRc6 


1.30 


210 


1.62x10* 


133 


2.69xl03 


IRc6 N 


0.87 


223 


2.09x10* 


133 


2.72x103 


IRc6 E 


1.48 


206 


1.52x10* 


134 


2.75 xl03 


IRc7 


2.02 


198 


1.14x10* 


148 


4.17x103 


IRc8 


0.79 


225 


2.16x10* 


119 


1.69x103 


IRcll 


1.73 


215 


1.87x10* 


144 


3.69x103 


IRcl2 


1.79 


258 


3.93x10* 


140 


3.34x103 


IRcl3 


0.84 


233 


2.49x10* 


118 


1.66x103 


IRcl4 


0.95 


223 


2.11x10* 


124 


2.03xl03 


IRcl5 


1.00 


222 


2.06x10* 


130 


2.41xl03 


IRcl6 N 


1.06 


238 


2.70x10* 


128 


2.32xl03 


IRcl6 S 


1.11 


227 


2.25x10* 


136 


2.97x10'^ 


IRcl6 W 


0.86 


233 


2.47x10* 


140 


3.38x103 


IRcl7 


0.65 


247 


3.16x10* 


122 


1.89xl03 


IRcig 


1.02 


303 


7.30x10* 


131 


2.48xl03 


IRc20 


1.33 


214 


1.75x10* 


126 


2.18x10-^ 


IRc21 


1.63 


221 


2.03x10* 


134 


2.74x10-^ 


IRc22 


1.05 


229 


2.36x10* 


126 


2.15x103 


IRc23 


1.21 


239 


2.80x10* 


126 


2.12xl03 


Knot 1 


0.81 


240 


2.83x10* 


115 


1.48xl03 


Knot 2 


0.84 


246 


3.12x10* 


114 


1.43x103 


Knot 3 


0.85 


245 


3.06x10* 


115 


1.50xl03 


BN SW 


0.77 


234 


2.53x10* 


140 


3.30 xl03 


MAX-69 




297 


6.79x10* 


128 


2.36 Xl03 



